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The  Modeling  of  Bistatic  Scattering  With  Moving  Platforms 


Abstract  ‘ 

The  complex  frequency  response  of  the  ocean  is  derived  for  three  different  bistatic  scattering 
problems.  The  derivations  are  based  on  treating  the  speed  of  sound  and  ambient  density  of  the 
ocean  as  constants,  and  solving  for  the  direct  ray  path  between  transmitter  and  discrete  point 
scatterer,  and  from  discrete  point  scatterer  to  receiver.  The  bistatic  scattering  problems  considered 
are:  1)  no  motion,  2)  only  the  discrete  point  scatterer  is  in  motion,  and  3)  all  three  platforms  (the 
transmitter,  discrete  point  scatterer,  and  receiver)  are  in  motion.  The  first  bistatic  scattering  problem 
yields  a  time-invariant,  space- variant  complex  frequency  response  while  the  remaining  two  bistatic 
scattering  problems  yield  time-variant,  space-variant  complex  frequency  responses.  For  problems 
involving  motion,  the  exact  time- varying  ranges  between  the  transmitter  and  discrete  point  scatterer, 
and  between  the  discrete  point  scatterer  and  receiver  are  derived,  and  the  exact  time-varying  angles 
of  incidence  at  the  discrete  point  scatterer,  and  the  exact  time-varying  angles  of  scatter  at  the  receiver 
are  also  derived.  The  solutions  for  the  exact  time- varying  ranges  are  also  valid  in  an  inhomogeneous 
ocean  where  the  speed  of  sound  and  ambient  density  are  functions  of  position  since  solving  for  a 
range  represents  a  problem  in  mechanics  not  wave  propagation. 


1  Introduction 

The  main  purpose  of  this  report  is  to  derive  the  complex  frequency  response  of  the  ocean  for 
three  different  bistatic  scattering  problems.  Since  we  will  be  working  with  small-amplitude  acoustic 
signals,  a  linear  wave  equation  accurately  describes  the  propagation  of  sound  between  source, 
discrete  point  scatterer,  and  receiver.  As  a  result,  we  treat  the  propagation  of  small-amplitude 
acoustic  signals  in  the  ocean  as  transmission  through  a  linear,  time-variant,  space-variant  filter. 
Treating  the  ocean  medium  as  a  linear  filter  is  valid  because  we  are  trying  to  solve  a  linear  wave 
equation. 

Scatter  from  a  discrete  point  scatterer  is  modeled  via  the  scattering  function,  which  is  a 
complex  function  (magnitude  and  phase)  and  is,  in  general,  a  function  of  frequency,  the  direction  of 
wave  propagation  from  the  source  to  the  scatterer,  and  the  direction  of  wave  propagation  from  the 
scatterer  to  the  receiver.  In  addition  to  the  scattering  function,  frequency-dependent  attenuation  is 
taken  into  account  in  order  to  model  the  propagation  of  sound  from  transmitter  to  discrete  point 
scatterer,  and  from  discrete  point  scatterer  to  receiver. 

The  speed  of  sound  and  ambient  density  of  the  ocean  are  treated  as  constants.  Therefore, 
sound  rays  will  travel  in  straight  lines.  We  will  only  concern  ourselves  with  solving  for  the  direct 
ray  path  between  transmitter  and  discrete  point  scatterer,  and  from  discrete  point  scatterer  to 
receiver.  As  a  result,  the  three  platforms  can  be  treated  as  being  in  an  unbounded,  homogeneous 
ocean  medium.  Sound  propagation  between  the  transmitter  and  the  ocean  surface  and  bottom,  and 
from  the  ocean  surface  and  bottom  to  the  receiver  can  be  handled  in  the  same  way  as  will  be 
developed  for  the  discrete  point  scatterer. 

Section  2.1  is  devoted  to  the  first  bistatic  scattering  problem,  which  involves  no  motion  -  the 
transmitter,  discrete  point  scatterer,  and  receiver  are  not  in  motion.  Most  of  Section  2.1  follows  the 
analysis  presented  in  [1].  However,  what  is  new  in  Section  2.1  is  the  derivation  of  the  exact  solution 
for  the  angles  of  incidence  at  the  discrete  point  scatterer  and  the  angles  of  scatter  at  the  receiver.  An 
example  is  worked  out  at  the  end  of  Section  2.1  showing  how  the  general  bistatic  scattering  results 
reduce  for  a  monostatic  ( backscatter )  scattering  geometry  for  the  no  motion  case. 

Section  2.2  is  devoted  to  the  second  bistatic  scattering  problem  when  only  the  discrete  point 
scatterer  is  in  motion.  The  initial  analysis  in  Section  2.2  follows  that  presented  in  [2],  with  the 
exception  that  motion  is  now  allowed  to  start  at  an  arbitrary  time  instant  tm  seconds  as  opposed  to 
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zero  seconds  as  was  done  in  [2].  In  addition,  two  new  major  results  are  presented  in  Section  2.2:  1) 
the  exact  time-varying  ranges  between  the  transmitter  and  discrete  point  scatterer,  and  between  the 
discrete  point  scatterer  and  receiver  are  derived,  and  2)  the  exact  time-varying  angles  of  incidence  at 
the  discrete  point  scatterer,  and  the  exact  time-varying  angles  of  scatter  at  the  receiver  are  also 
derived.  An  example  is  worked  out  at  the  end  of  Section  2.2  showing  how  the  general  bistatic 
scattering  results  reduce  for  a  monostatic  (backscatter)  scattering  geometry  for  the  case  when  only 
the  discrete  point  scatterer  is  in  motion. 

Section  2.3  is  devoted  to  the  third  bistatic  scattering  problem  when  all  three  platforms  are  in 
motion.  The  initial  analysis  in  Section  2.3  follows  that  presented  in  [3],  with  the  exception  that 
motion  is  now  allowed  to  start  at  an  arbitrary  time  instant  tm  seconds  as  opposed  to  zero  seconds  as 
was  done  in  [3].  In  addition,  two  new  major  results  are  presented  in  Section  2.3:  1)  the  exact  time- 
varying  ranges  between  the  transmitter  and  discrete  point  scatterer,  and  between  the  discrete  point 
scatterer  and  receiver  are  derived,  and  2)  the  exact  time-varying  angles  of  incidence  at  the  discrete 
point  scatterer,  and  the  exact  time-varying  angles  of  scatter  at  the  receiver  are  also  derived.  Three 
examples  are  worked  out  at  the  end  of  Section  2.3.  The  first  example  shows  that  the  exact  results 
derived  in  Section  2.3  reduce  to  the  exact  results  derived  in  Sections  2.1  and  2.2  when  appropriate 
values  are  used  for  the  various  parameters.  This  is  a  very  important  example  because  it  validates  the 
correctness  of  the  general  solution  derived  in  Section  2.3.  The  second  example  shows  how  the 
general  bistatic  scattering  results  reduce  for  a  monostatic  ( backscatter )  scattering  geometry  for  the 
case  when  all  three  platforms  are  in  motion.  The  third  example  shows  how  die  general  bistatic 
scattering  results  can  be  applied  to  a  synthetic  aperture  sonar  (SAS)  system  trying  to  image  a 
nonmoving  target  on  the  ocean  bottom  without  having  to  make  several  simplifying  assumptions  as 
is  done,  for  example,  in  [8]. 

Finally,  note  that  for  problems  involving  motion,  the  solutions  for  the  exact  time-varying 
ranges  between  the  transmitter  and  discrete  point  scatterer,  and  between  the  discrete  point  scatterer 
and  receiver  derived  in  this  report  are  also  valid  in  an  inhomogeneous  ocean  where  the  speed  of 
sound  and  ambient  density  are  functions  of  position  since  solving  for  a  range  represents  a  problem 
in  mechanics  not  wave  propagation.  However,  travel  times  and  angles  of  incidence  and  scatter  are 
different  in  an  inhomogeneous  ocean  compared  to  a  homogeneous  ocean  because  of  the 
complicated  trajectories  of  sound  rays  in  an  inhomogeneous  ocean. 


2  Bistatic  Scattering 
2.1  No  Motion 

In  this  section  we  will  analyze  the  simple  bistatic  scattering  problem  shown  in  Fig.  2.1-1.  The 
transmitter,  discrete  point  scatterer,  and  receiver  are  not  in  motion.  No  motion  corresponds  to  a 
time-invariant  problem.  As  mentioned  in  the  Introduction,  the  three  platforms  will  be  treated  as 
being  in  an  unbounded,  homogeneous  ocean  medium.  Although  the  propagation  of  sound  between 
the  source  and  discrete  point  scatterer,  and  between  the  discrete  point  scatterer  and  receiver  can  be 
treated  as  transmission  through  linear,  time-invariant,  space-invariant  filters;  the  overall  solution 
for  this  bistatic  scattering  problem  corresponds  to  transmission  through  a  linear,  time-invariant, 
space-variant  filter.  The  presence  of  a  discrete  point  scatterer  in  an  unbounded,  homogeneous  fluid 
medium  (i.e.,  a  fluid  medium  with  constant  speed  of  sound  and  ambient  density)  causes  the  medium 
to  be  space- variant. 

Let  the  source  distribution  xM(t,  r)  at  time  t  and  position  r  =  (x,y,z)  be  a  motionless,  time- 
harmonic,  point  source  with  units  of  inverse  seconds,  that  is,  let 

xM(t,r)  =  S0S(r-r0)exp(+j2nft),  (2.1-1) 
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where  S0  is  the  source  strength  in  cubic  meters  per  second,  the  impulse  function  <5(r-r0),  with 
units  of  inverse  cubic  meters,  represents  a  point  source  at  r0  =  (x0,y0,z0)’ anc*  /  is  frequency  in 
hertz.  The  sound  source  has  been  turned  on  forever,  that  is,  since  t  =  -°°.  The  propagation  of 
sound  between  the  source  and  discrete  point  scatterer  can  be  modeled  as  transmission  through  a 
linear,  time-invariant,  space -invariant,  filter.  Therefore,  the  acoustic  field  (velocity  potential) 
incident  upon  the  discrete  point  scatterer  at  r,  =  (x^y,,*,)  is  given  by  [4] 

yM(t,  r,)  =  yfM(ri)exp(+j2nft),  (2.1-2) 

where 


y ^  m  ( y*  ^*i  r0). 


H, 


M 


,,  x  exp(-^r,-r0) 
(/,r,-r0)  = - 

r,-r0 


(2.1-3) 

(2.1-4) 


is  the  time-invariant,  space-invariant,  complex  frequency  response  of  the  ocean  at  frequency  / 
hertz, 


k  -  2nf  /c  =  2n/X  (2.1-5) 

is  the  wavenumber  with  units  of  radians  per  meter,  c  is  the  constant  speed  of  sound  of  the 
homogeneous  ocean  medium  in  meters  per  second,  and  A  is  the  wavelength  in  meters. 

As  can  be  seen  from  Fig.  2.1-1,  the  position  vector  from  the  point  source  to  the  discrete  point 
scatterer  is  given  by 


Therefore,  (2.1-3)  and  (2.1-4)  can  be  rewritten  as 


—  ‘■*0  ( •/"’  ®"o,i ) 


(2.1-6) 


(2.1-7) 


and 


HM(f,r0A)  =  - 


exp| 

Ml' 

r»'l) 

An 

r0. 

,1 

(2.1-8) 


respectively. 

In  order  to  compute  the  acoustic  signal  incident  upon  the  receiver,  we  treat  the  discrete  point 
scatterer  as  another  motionless,  time-harmonic,  point  source  with  units  of  inverse  seconds,  that  is, 
let  [see  (2.1-1)  and  Fig.  2.1-2] 


x'M(t,  r)  =  S'<5(r  -  r,)exp(+;2^/0,  (2.1-9) 

where  S'Q  is  the  source  strength  in  cubic  meters  per  second,  and  the  impulse  function  <5(r-r,), 
with  units  of  inverse  cubic  meters,  represents  a  point  source  at  r,  =  (x,,y,,z,).  The  source  strength 
S'}  is  given  by 
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(2.1-10) 


where  y/iM(r,)  is  given  by  (2.1-3)  and  is  the  spatial-dependent  part  of  the  time-harmonic  velocity 
potential  at  r,  with  units  of  squared-meters  per  second,  f,n0A,nl2)  is  the  scattering  function  of 
the  discrete  point  scatterer  with  units  of  meters,  and  n0 ,  and  hl2  are  the  dimensionless  unit  vectors 
in  the  directions  measured  from  the  source  to  the  discrete  point  scatterer,  and  from  the  discrete  point 
scatterer  to  the  receiver,  respectively.  The  unit  of  meters  for  the  scattering  function  represents  an 
effective  scattering  length  that  may  be  larger  or  smaller  than  the  actual  length  of  the  scatterer.  The 
scattering  function  is  a  complex  function  (magnitude  and  phase)  and  is,  in  general,  a  function  of 
frequency,  the  direction  of  wave  propagation  from  the  source  to  the  scatterer  ( n01),  and  the  direction 
of  wave  propagation  from  the  scatterer  to  the  receiver  (n,  2)  (see  either  Fig.  2.1-1  or  Fig.  2.1-2)  [5]. 
Later  in  this  section,  we  will  show  how  to  express  the  scattering  function  as  not  only  a  function  of 
frequency,  but  also  as  a  function  of  the  angles  of  incidence  and  scatter  instead  of  unit  vectors.  The 
use  of  unit  vectors  is  meant  as  a  shorthand  notation. 

The  propagation  of  sound  between  the  discrete  point  scatterer  and  receiver  can  also  be 
modeled  as  transmission  through  a  linear,  time-invariant,  space-invariant  filter.  Therefore,  the 
acoustic  field  (velocity  potential)  incident  upon  the  receiver  at  r2  =  (x2,y2,z2)  due  to  a  point  source 
at  r,  =  (x,,y,,z,)  is  given  by  [4] 


yM(t,r2)  =  yfM(r2)exp(+j2nft), 


where 


and 


7/.*  (F2  )  =  ^o' ***(/,  Ij-r,) 


Hu(f,r2-  r,) 


exp(-,/*|r2-r1|) 

47r|r2-r,| 


(2.1-11) 


(2.1-12) 


(2.1-13) 


As  can  be  seen  from  either  Fig.  2.1-1  or  Fig.  2.1-2,  the  position  vector  from  the  discrete  point 
scatterer  to  the  receiver  is  given  by 


Therefore,  (2.1-12)  and  (2.1-13)  can  be  rewritten  as 


y/,  M  (r2  )  —  S0  ( /,  r,  2 ) 


(2.1-14) 


(2.1-15) 


and 


HM(f,rl2)  =  - 


exp| 

b'4 

ru|) 

An 

ri.: 

>  1 

respectively. 


(2.1-16) 
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Let  us  now  begin  the  process  of  obtaining  a  final  expression  for  the  time-harmonic  velocity 
potential  incident  upon  the  receiver.  Substituting  (2.1-10)  and  (2.1-7)  into  (2.1-15)  yields 


or,  equivalently, 


where 


}'/.M(r2)  =  ^0^Af(/>r0.l)Sl(/>"o.P«1.2)Ww(/»ri.2)’ 

y /,m  (®2 ) =  S0HM(f,r2 1  r0), 

Hm ( /.I r2  |  ro )  =  Hm ( /. : 1*0.1 ) 8 1  ( fj «0,P «1.2 ) #  m ( /. **1.2 ) 

< 

=  «i(/.«o.iA2)- 
=  8i{f>n o.i^u)' 


exp 

~jk( 

ro,i 

+ 

M)] 

16tt2 

ro,i 

ri,: 

>| 

exp 

~jk( 

ri_ro 

+  r2  “  ri  ) 

16^2 

ri-ro 

r2-r, 

(2.1-17) 


(2.1-18) 


(2.1-19) 


Note  that  if  the  bistatic  scattering  problem  shown  in  Fig.  2.1-1  corresponded  to  transmission 
through  a  space-invariant  filter,  then  the  complex  frequency  response  given  by  (2.1-19)  would  be  a 
function  of  the  vector  spatial  difference  r2  -  r0,  which  it  is  not.  Also  note  that  in  order  to  model  the 
effects  of  frequency-dependent  attenuation,  simply  replace  the  real  wavenumber  k  given  by  (2.1-5) 
with  the  following  complex  wavenumber  K : 


K  =  k-ja(f), 


(2.1-20) 


where  cc(f)  is  the  real,  frequency-dependent,  attenuation  coefficient  in  nepers  per  meter.  In 
addition  to  being  real  quantities,  both  k  and  a(f)  are  positive. 

With  the  use  of  (2.1-5),  (2.1-11),  (2.1-18),  and  (2.1-19);  and  by  replacing  the  real 
wavenumber  k  in  (2.1-19)  with  the  complex  wavenumber  K  given  by  (2.1-20),  we  can  summarize 
our  results  as  follows:  for  the  bistatic  scattering  problem  shown  in  Fig.  2.1-1,  the  time-harmonic 
velocity  potential  in  squared-meters  per  second  incident  upon  the  receiver  at  r2  =  (x2,y2,z2),  due  to 
a  time-harmonic  point  source  at  r0  =  (x0,y0,z0)  and  a  discrete  point  scatterer  at  r,  =  (Xpy^Zj),  is 
given  by 


yM(t,r2)  =  S0Hm(  f,  r2 1  r0)exp(+;2^/O,  t  >  r, 


(2.1-21) 


where  S0  is  the  source  strength  in  cubic  meters  per  second, 


M  (  f  v  \r\—  o  (  f  ri  fi  ^ 

-«(/)('■«, +n.1)] 

/’r2  I  ro )  ~  S\\J / 

2  eXP^  JLJLJ  L ) 

r01r12 

(2.1-22) 
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is  the  time-invariant,  space-variant,  complex  frequency  response  of  the  ocean  at  frequency  /  hertz, 
5 1(  M.1A2)  is  scattering  function  of  the  discrete  point  scatterer  in  meters  [see  (2.1-46)], 


and 


(2.1-23) 


(2.1-24) 


are  the  dimensionless  unit  vectors  in  the  directions  measured  from  the  source  to  the  discrete  point 
scatterer,  and  from  the  discrete  point  scatterer  to  the  receiver,  respectively,  a(f)  is  the  real, 
frequency-dependent,  attenuation  coefficient  in  nepers  per  meter, 


(2.1-25) 


is  the  time  delay  in  seconds  (the  amount  of  time  it  takes  for  the  transmitted  acoustic  signal  to  begin 
to  appear  at  the  receiver),  and 


ro,i  = 


(2.1-26) 


and 


(2.1-27) 


are  the  distances  (ranges)  in  meters  measured  from  the  source  to  the  discrete  point  scatterer,  and 
from  the  discrete  point  scatterer  to  the  receiver,  respectively. 

Equations  (2.1-23)  and  (2.1-24)  represent  one  way  to  compute  the  unit  vectors  n0I  and  hl2. 
Unit  vectors  can  also  be  computed  using  direction  cosines.  For  example,  the  dimensionless  unit 
vector  n01  can  also  be  expressed  as  follows: 


where 


and 


«0,l^+V0.li)  +  W0.1^ 

(2.1-28) 

i  —  sin  Oq  j  cos  i jf  Qti  9 

(2.1-29) 

,  =  sin0ol  siny'oj, 

(2.1-30) 

W0l  =  COS  0Q.1 

(2.1-31) 
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are  dimensionless  direction  cosines  with  respect  to  the  X,  Y,  and  Z  axes,  respectively,  and  0O1  and 
y/0i  are  spherical  angles  defined  in  Fig.  2.1-3.  Equating  (2.1-23)  and  (2.1-28)  yields 


and 


M0,l~ 


*i-*o 

r,-r0 


(2.1-32) 


V0,l  = 


yi-y0 

ki-r0 


(2.1-33) 


(2.1-34) 

lri-ro| 

Therefore,  with  the  use  of  (2.1-29)  through  (2.1-31),  and  (2.1-32)  through  (2.1-34),  the  spherical 
angles  0O ,  and  y/0  l  can  be  computed  as  follows: 


0ol  =  cos  1  w01  =  cos 


(2.1-35) 


and 


(2.1-36) 


Equations  (2.1-35)  and  (2.1-36)  are  the  angles  of  incidence  at  the  discrete  point  scatterer. 
Similarly,  the  dimensionless  unit  vector  hl2  can  also  be  expressed  as  follows: 


\2=ui,2^ +  vl2y  +  wl2z, 

where 

1^2=  sin  0{  2  cos  y/l  2 , 
v,  2=  sin  012  sin^12» , 


and 


(2.1-37) 


(2.1-38) 

(2.1-39) 


W,  2=  COS  0,  2  (2.1-40) 

are  dimensionless  direction  cosines  with  respect  to  the  X,  Y,  and  Z  axes,  respectively,  and  0l2  and 
y/l  2  are  spherical  angles  defined  in  Fig.  2.1-4.  Equating  (2.1-24)  and  (2.1-37)  yields 
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z 


Figure  2.1-3.  The  dimensionless  unit  vector  n0l  measured  in  the  direction  from  the  source  to  the 
discrete  point  scatterer,  and  associated  spherical  angles  60l  and  \j/0l. 
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z 


Figure  2.1-4.  The  dimensionless  unit  vector  hl  2  measured  in  the  direction  from  the  discrete  point 
scatterer  to  the  receiver,  and  associated  spherical  angles  d]  2  and  \j/l2. 
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(2.1-43) 


Therefore,  with  the  use  of  (2.1-38)  through  (2.1-40),  and  (2.1-41)  through  (2.1-43),  the  spherical 
angles  0, 2  and  \ffi  2  can  be  computed  as  follows: 


and 


(2.1-45) 


Equations  (2.1-44)  and  (2.1-45)  are  the  angles  of  scatter  at  the  receiver. 

The  result  of  the  above  analysis  is  that  the  scattering  function  can  also  be  expressed  as  a 
function  of  two  sets  of  angles  -  the  angles  of  incidence  (0OJ,v^o.i)  given  by  (2.1-35)  and  (2.1-36), 
and  the  angles  of  scatter  (012,i/r12)  given  by  (2.1-44)  and  (2.1-45).  Therefore, 

(2.1-46) 

Let  us  next  relate  the  scattering  function,  the  differential  scattering  cross  section,  and  target 
strength  of  the  discrete  point  scatterer.  The  target  strength  (TS)  is  defined  as  follows  [6]: 


(2.1-47) 


where  [5, 6] 
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(2.1-48) 


®d{  /’^ 0,1 » 


ri,2  Kvgs  (*2) 

_  S.(/Ap”u)| 

LgM 1) 

(47T)2 

is  the  differential  scattering  cross  section  with  units  of  squared  meters,  I  (r,)and  7  ( r2)  are  the 

time-average,  incident  and  scattered  intensities,  respectively,  with  units  of  watts  per  squared  meter, 
£i(/>"o,p”i,2)  is  ^e  scattering  function  of  the  discrete  point  scatterer  with  units  of  meters,  and  Aref 
is  a  reference  cross-sectional  area  commonly  chosen  to  be  equal  to  1  m2 . 

Example  2.1-1  Monostatic  Scattering  Geometry 

For  a  monostatic  ( backscatter )  scattering  geometry,  both  the  transmitter  and  receiver  are 
located  at  the  same  position,  that  is, 

r2  =  r0.  (2.1-49) 

Substituting  (2.1-49)  into  (2.1-14)  yields 

ri.2  =  “ro,P  (2.1-50) 

where  r01  is  given  by  (2.1-6).  Therefore, 

r12  =  r0,  .  (2.1-51) 

and 

«,2  =  -n0,.  (2.1-52) 

With  the  use  of  (2.1-51)  and  (2.1-52),  the  time-invariant,  space-variant,  complex  frequency 
response  of  the  ocean  given  by  (2.1-22)  reduces  to 


where 


Tj  t  r  i  \  /,  *  -  \  exp[-2a(/)r01] 

hm  ( M I  r0 )  =  g ,  ( /,  n0>1 ,  ■ -  n0,i ) - 77 - 75 - -  exp  (-j  2nff) , 

(47Tr0  j ; 


(2.1-53) 


^0,1)  —  £i( + 


(2.1-54) 


(60l,y/ol)  are  the  angles  of  incidence  given  by  (2.1-35)  and  (2.1-36), 
V\,2  +  Vo.i  31:6  ^  angles  of  scatter. 


012  =  n  -  Q0 ,  and 
(2.1-55) 


is  the  time  delay  in  seconds,  and  r0 ,  is  given  by  (2.1-26). 
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2.2  Discrete  Point  Scatterer  In  Motion 

In  this  section  we  will  analyze  the  bistatic  scattering  problem  shown  in  Fig.  2.2-1.  The 
transmitter  and  receiver  are  not  in  motion.  Only  the  discrete  point  scatterer  is  in  motion.  Motion 
corresponds  to  a  time-variant  problem.  As  mentioned  in  the  Introduction,  the  three  platforms  will 
be  treated  as  being  in  an  unbounded,  homogeneous  ocean  medium.  Although  the  propagation  of 
sound  between  the  source  and  discrete  point  scatterer,  and  between  the  discrete  point  scatterer  and 
receiver  can  be  treated  as  transmission  through  linear,  time-variant ,  space-invariant  filters;  the 
overall  solution  for  this  bistatic  scattering  problem  corresponds  to  transmission  through  a  linear, 
time-variant,  space-variant  filter.  The  presence  of  a  discrete  point  scatterer  in  an  unbounded, 
homogeneous  fluid  medium  (i.e.,  a  fluid  medium  with  constant  speed  of  sound  and  ambient 
density)  causes  the  medium  to  be  space-variant. 

Let  the  source  distribution  xM(t, r)  at  time  t  and  position  r  =  (jc,y,z)  be  a  motionless,  time- 
harmonic,  point  source  with  units  of  inverse  seconds,  that  is,  let 

xM(t,  r)  =  S0<5(r-ro)exp(+;27r/O,  (2.2-1) 

where  S0  is  the  source  strength  in  cubic  meters  per  second,  the  impulse  function  <5(r  -  r0),  with 
units  of  inverse  cubic  meters,  represents  a  point  source  at  r0  =  (jr0,y0,z0),  and  /  is  frequency  in 
hertz.  The  sound  source  has  been  turned  on  forever,  that  is,  since  t  =  —  °°.  The  velocity  vector  of  the 
discrete  point  scatterer,  V,,  is  given  by 


V,=V;«v  (2.2-2) 

where  V]  is  the  speed  in  meters  per  second  and  hv  is  the  dimensionless  unit  vector  in  the  direction 
of  V, .  The  velocity  vector  given  by  (2.2-2)  is  constant,  that  is,  both  the  speed  and  direction  are 
constants  -  there  is  no  acceleration.  Motion  begins  at  time  t—tm  seconds.  We  will  model  the 
propagation  of  sound  from  the  time  motion  begins. 

Since  the  discrete  point  scatterer  is  now  in  motion,  the  position  vector  from  the  origin  to  the 
discrete  point  scatterer,  denoted  by  )?,(/),  is  a  function  of  time  given  by 

»I(0  =  r,  +  A/V„  t>tm,  (2.2-3) 

where  r,  =(x,,yl,zI)  is  the  position  vector  from  the  origin  to  the  discrete  point  scatterer  when 
motion  begins  (see  Fig.  2.2-1),  and 


A  t  =  t-t  t> t 


(2.2-4) 


Note  that  /^y(tm)~  r,. 

When  the  transmitted  acoustic  field  is  first  incident  upon  the  discrete  point  scatterer  at  some 
time  t'  seconds  where  t'  >  tm,  the  position  vector  from  the  origin  to  the  discrete  point  scatterer  is 
given  by  [see  (2.2-3)  and  Fig.  2.2-2] 

r,'  =£,i(f')  =  r1  +  A/'  V, ,  t'>tm,  (2.2-5) 

where 

A t'=t'-tm,  t'  >tm.  (2.2-6) 
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Figure  2.2-1.  Bistatic  scattering  geometry  when  motion  begins  at  time  t=tm  seconds.  Point  0, 
P0(r0),  is  the  transmitter;  point  1,  /{(r,),  is  the  discrete  point  scatterer;  and  point  2,  P2(r2),  is  the 
receiver.  Only  the  discrete  point  scatterer  is  in  motion. 
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««(!•- O 


S0S(  r- 


«1.2 


Figure  2.2-2.  Bistatic  scattering  geometry  when  the  transmitted  acoustic  field  is  first  incident 
upon  the  discrete  point  scatterer  at  time  t’  seconds.  Point  0,  P0(r0),  is  the  transmitter;  point  1, 
/{(r,),  is  the  discrete  point  scatterer;  and  point  2,  P2(r2),  is  the  receiver.  Only  the  discrete  point 
scatterer  is  in  motion. 
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The  propagation  of  sound  between  the  source  and  discrete  point  scatterer  can  be  modeled  as 
transmission  through  a  linear,  time-variant,  space -invariant  filter.  Therefore,  the  acoustic  field 
(velocity  potential)  incident  upon  the  discrete  point  scatterer  at  time  t'  and  position  r,'  =  {x[,y[,z[ ) 
is  given  by  [7] 


yM  (/',  <  )  =  S0Hm  (t\ r(  -  r0 1  /)  exp(+j2nft' ) ,  t'>tm,  (2.2-7) 


where 


exp(-;fc|r;-r0|) 

4^|r;-r0| 


(2.2-8) 


is  the  time-variant,  space-invariant,  complex  frequency  response  of  the  ocean  at  frequency  /  hertz, 
and  k  is  the  wavenumber  in  radians  per  meter  given  by  (2.1-5)  and  is  repeated  below  for 
convenience: 


k  =  lizfjc  =  In/X. 


(2.1-5) 


By  referring  to  Fig.  2.2-2,  we  can  express  the  position  vector  from  the  point  source  to  the 
discrete  point  scatterer  at  time  t'  as 


ro.i  =  ri  ~  ro> 

and  upon  substituting  (2.2-5)  into  (2.2-9),  we  obtain 

r„,i =  r0Ai =  **o,i  ^ 


(2.2-9) 


(2.2-10) 


where  r'0l  = 


0,1 


and  is  repeated 


,  h'0 ,  is  the  dimensionless  unit  vector  in  the  direction  of 
jelow  for  convenience  (also  see  Fig.  2.2-2), 


L0,P 


r0  l  is  given  by  (2.1-6) 


ro,i=  ri-I0’ 


(2.1-6) 


and  At'  is  given  by  (2.2-6).  Therefore,  (2.2-7)  and  (2.2-8)  can  be  rewritten  as 
y m )  =  S0Hm  (t' ,r',| /) exp(+j2nft'),  t’>tm. 


and 


(2.2-11) 


(2.2-12) 


In  order  to  compute  the  acoustic  signal  incident  upon  the  receiver,  we  treat  the  discrete  point 
scatterer  at  time  t>t'>tm  and  position  r'  as  another  motionless,  time-harmonic,  point  source  with 
units  of  inverse  seconds,  that  is,  let  [see  (2.2-1)  and  Fig.  2.2-2] 

x'M(t,r)  =  S'08(r-r')exp(+j2nft),  (2.2-13) 
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where  Sq  is  the  source  strength  in  cubic  meters  per  second  and  the  impulse  function  <5(r  -  r,' ), 
with  units  of  inverse  cubic  meters,  represents  a  point  source  at  r,'  =  (x[,y[,z[ ) .  The  source  strength 
Sq  is  given  by 


(2.2-14) 


where  g ,  ( /,  h'0 , ,  h[  2 )  is  the  scattering  junction  of  the  discrete  point  scatterer  with  units  of  meters  as 
discussed  in  Section  2.1.  Since  the  discrete  point  scatterer  is  now  in  motion  and  the  transmitted 
acoustic  field  is  first  incident  upon  the  discrete  point  scatterer  at  time  t',  the  direction  of  wave 
propagation  from  the  source  to  the  scatterer  is  given  by  the  dimensionless  unit  vector  h'0 ,  (see  Fig. 
2.2-2).  Similarly,  since  the  discrete  point  scatterer  is  being  treated  as  another  point  source  at  time  t' 
and  position  r,' ,  the  direction  of  wave  propagation  from  the  scatterer  to  the  receiver  is  given  by  the 
dimensionless  unit  vector  h[2  (see  Fig.  2.2-2).  Later  in  this  section,  we  will  show  that  the  scattering 
function  is  also  a  function  of  time  because  the  unit  vectors  are  actually  time-varying.  We  will  then 
show  how  to  express  the  scattering  function  as  a  function  of  frequency  and  time-varying  angles  of 
incidence  and  scatter  instead  of  time-varying  unit  vectors.  The  use  of  unit  vectors  is  meant  as  a 
shorthand  notation. 

The  propagation  of  sound  between  the  discrete  point  scatterer  and  receiver  can  also  be 
modeled  as  transmission  through  a  linear,  time-variant,  space-invariant  filter.  The  scattered 
acoustic  field  is  first  incident  upon  the  receiver  at  some  time  t  seconds  where  t>  t'  >  tm.  Therefore, 
the  acoustic  field  (velocity  potential)  incident  upon  the  receiver  at  time  t  and  position  r2  =  (x2,  y2,z2) 
due  to  a  point  source  at  r,'  =  (*,',  y[,  z[ )  is  given  by  [7] 


where 


yM(t,r2)=  S'0HM(t,r2-r;\f)exp(+j27ift),  t>t'>tm,  (2.2-15) 


exp(-^|r2-<|) 

4rc|r2-r1'| 


(2.2-16) 


By  referring  to  Fig.  2.2-2,  we  can  express  the  position  vector  from  the  discrete  point  scatterer 
to  the  receiver  at  time  t'  >  tm  as 


—  r  —  r/ 

1 1,2  l2  *1  ® 

and  upon  substituting  (2.2-5)  into  (2.2-17),  we  obtain 

'  '  _  A  t'\ 

*  1,2  '1,2  '*1,2  *  1,2  m  VP 


(2.2-17) 


(2.2-18) 


where  r'2=  r,' 


1 1.2 


,  h'i  2  is  the  dimensionless  unit  vector  in  the  direction  of  r'2,  r12  is  given  by  (2.1- 


14)  and  is  repeated  below  for  convenience  (also  see  Fig.  2.2-2), 

r„ 


ri,2  =  r2  ‘ 


(2.1-14) 


and  A t'  is  given  by  (2.2-6).  Therefore,  (2.2-15)  and  (2.2-16)  can  be  rewritten  as 
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y  m  (t,r2)  =  S'ohm  (t,r'l2\f)exp(+j2nft),  t>t'>tm,  (2.2-19) 

and 


H, 


exp| 

M1 

ru|) 

An 

r[, 

1*^ 

)  1 

(2.2-20) 


Let  us  now  begin  the  process  of  obtaining  a  final  expression  for  the  time-harmonic  velocity 
potential  incident  upon  the  receiver.  Substituting  (2.2-14)  into  (2.2-19)  yields 


y  m  (^»r2) —  S0Hm  1 t  ,r0,|/jg,(  f,nQvnh2)HM  |t,r12|/jexp(-t-_/2tr/t), 


t>t'  >t. 


m’ 


(2.2-21) 


or,  equivalently, 

y^n (^t*2) =  SqH M (t,r2  |y",r0)exp(+y'2?Ti/'f),  t>t  >  tm,  (2.2-22) 

where 


—  &l(  /’^O.l’^u) 
=  £i(  /»«0,P«1,2) 


exp 

~jk[ 

*0.1 

+ 

<>l)] 

167T2 

r0,i 

K: 

J 

exp[-./fc(|r1  -  r0  +  At' Vj |  +  |r2  -  r,  -  At' V,|) 
167r2|r,  -  r0  +  At'V;  |r2  -  ^  -  At'  V,| 


At'  is  given  by  (2.2-6),  and 


t  =  t'  +  ■ 


1 1,2 1 


(2.2-23) 


(2.2-24) 


or 


(2.2-25) 
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where  r'l2  is  given  by  (2.2-18).  Note  that  if  the  bistatic  scattering  problem  shown  in  Fig.  2.2-1 
corresponded  to  transmission  through  a  space-invariant  filter,  then  the  complex  frequency  response 
given  by  (2.2-23)  would  be  a  function  of  the  vector  spatial  difference  r2  -  r0,  which  it  is  not. 

In  order  to  evaluate  the  complex  frequency  response  given  by  (2.2-23),  we  must  derive 


solutions  for  the  unit  vectors  n 
quantities.  Let  us  begin  with 


and  n',2,  and  the  ranges 


0,1 


0.1 

.  Since 


0,1 


and  r 


‘1,2 


in  terms  of  known 


•  r, 


0.1  ’ 


(2.2-26) 


substituting  (2.2-10)  into  (2.2-26)  yields 

k,f  =  (r0>1  +  Af' V.)  •  (r01  +  At' V,).  (2.2-27) 


Expanding  the  right-hand  side  of  (2.2-27)  and  taking  the  square  root  of  both  sides  of  the  resulting 
equation  yields 


/  \ 

ro 
_ 1 

1 

1 

2 

Vl  At 

|r0,i 

II 

1  +  — (n01«  V.)Af  + 

{  'I,.,  J 

(2.2-28) 


or 


. 

V 

< 

i 

x  1  V,At' 

k. 

=  ro.f 

1  +  2— - 

(n0A»nv)  + 

ro.i 

2  r0A 

(2.2-29) 


Although  (2.2-28)  and  (2.2-29)  are  exact  expressions  for  the  range  r^,  ,  we  cannot  compute 

a  value  for  ,|  until  we  derive  an  expression  for  At'  in  terms  of  ,|.  By  referring  to  Fig.  2.2-2,  it 
can  be  seen  that 


t  >  t 


(2.2-30) 


Substituting  (2.2-30)  into  (2.2-28)  and  squaring  both  sides  of  the  resulting  equation  yields  the 
following  second-order  polynomial 


A 


-C  =  0, 


(2.2-31) 


with  exact  solution 
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(2.2-32) 


B±J¥+4AC 

_ 


where 


and 


A  =  1- 


f  Vj^2 


B  =  2r, 


\  c. 

VV, 


0,1 


(2.2-33) 

(2.2-34) 


C=  r2 
u  'o,i  • 


(2.2-35) 


The  solution  given  by  (2.2-32)  is  the  constant  value  of  range  between  the  point  source  and  the 
discrete  point  scatterer  when  the  transmitted  acoustic  field  is  first  incident  upon  the  discrete  point 
scatterer  at  time  instant  t'  after  motion  begins  at  time  instant  tm  where  t'  >  tm.  The  decision  to  use 
either  the  plus  or  minus  sign  in  (2.2-32)  is  dictated  by  the  fact  that  range  must  be  positive.  Let  us 

solve  for  |r,'2|  next. 

Since 


•  r; 


1,2* 


(2.2-36) 


substituting  (2.2-18)  into  (2.2-36)  yields 

|r'2|2  =  (r12  -  At'  V,)  •  (r12  -  At'  V,).  (2.2-37) 

Expanding  the  right-hand  side  of  (2.2-37)  and  taking  the  square  root  of  both  sides  of  the  resulting 
equation  yields 


or 


f 

2“ 

Ki\=hi 

1-—  (n12»V.)Ar,+ 

VxAt' 

1 

ro 

l  ru  J 

(2.2-38) 


=  r, 


1,2 


hi 


(nl2*nVi)~ 


1  VxAt' 

2  hi 


(2.2-39) 
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Although  (2.2-38)  and  (2.2-39)  are  exact  expressions  for  the  range  jr(2| ,  we  cannot  compute  a 
value  for  r[2  until  we  derive  an  appropriate  expression  for  A  t' .  We  will  first  solve  for  a  constant 
value  for  r'2  by  expressing  At'  in  terms  of  the  known  constant  range  |r„  ,|  [see  (2.2-32)].  We  will 
then  solve  for  a  time-varying  r,'2|  by  expressing  At' exclusively  in  terms  of  r[2|. 

Since  (2.2-30)  already  expresses  At'  in  terms  of  |r',|,  substituting  (2.2-30)  into  (2.2-38) 
yields 


(2.2-40) 


where  |r^  ,J  is  given  by  (2.2-32).  The  solution  given  by  (2.2-40)  is  the  constant  value  of  range 

between  the  discrete  point  scatterer  and  the  receiver  when  the  scattered  acoustic  field  is  first  incident 
upon  the  receiver  at  time  instant  t  after  motion  begins  at  time  instant  tm  where  t>  t'  >  tm. 

As  we  previously  mentioned,  we  will  now  solve  for  a  time- varying  |r,'2  by  expressing  At' 
exclusively  in  terms  of  r,'2|.  We  begin  by  substituting  (2.2-25)  into  (2.2-6)  which  yields 

|* 9  I 

A  t'=t-tm-^,  t>t'>tm,  (2.2-41) 

c 

and  upon  substituting  (2.2-4)  into  (2.2-41),  we  obtain 

(2.2-42) 

The  solution  for  |r,'2|  can  be  obtained  by  substituting  (2.2-42)  into  (2.2-38)  and  squaring 
both  sides  of  the  resulting  equation.  Doing  so  yields  the  following  second-order  polynomial 

^|r'2(/)|2  +  S(t)\r'l2(t)\- m  =  0,  t>tm  +  r,  (2.2-43) 

with  exact  solution 

(2.2-44) 

where 
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,  fvY 

^  =  i-  , 

(2.2-45) 

\  c  J 

V2  n,,*V, 

s 

ts 

sc 

CN 

1 

<1 

<N 

II 

/ — s 

•♦-a 

to 

(2.2-46) 

c  c 

m  =  V2(At)2-2rl2(nl2»  V,)A t+r?A, 

(2.2-47) 

A  t=t-tm,  t>tm  +  T, 

(2.2-48) 

and  (see  Fig.  2.2-2) 

(2.2-49) 

is  the  time  delay  in  seconds  (the  amount  of  time  it  takes  for  the  transmitted  acoustic  signal  to  begin 
to  appear  at  the  receiver  after  motion  begins  at  time  instant  tm )  where  the  constant  values  of  range 

[ Tq  ,  |  and  r,'2  are  given  by  (2.2-32)  and  (2.2-40),  respectively.  It  is  important  to  note  that  if  (2.2- 
43)  is  evaluated  at  t  =  tm  +  r,  then  it  can  be  shown  that 

<A  +  T)|  =  |r'2|,  (2.2-50) 

where  |r'2|  is  given  by  (2.2-40). 

Now  that  we  have  an  exact  solution  for  |r1'2(?)|  as  given  by  (2.2-44),  we  can  use  it  to  obtain  an 

exact  solution  for  the  time- varying  range  |r£j(f)|  as  follows.  With  the  use  of  (2.2-30)  and  (2.2-42), 
we  can  write  that 

|r;,(/)|  =  cAt-|r'2(0|,  t>tm  +  z,  (2.2-51) 

where  At  is  given  by  (2.2-48),  |r1'2(o|  is  given  by  (2.2-44),  and  r  is  given  by  (2.2-49).  Equation 
(2.2-51)  indicates  that  if  we  are  given  a  value  of  At,  which  determines  the  value  of  |r,'2(t)|,  then  we 
can  use  those  two  values  to  compute  what  r  (,,(/)  must  have  been  -  we  are  working  backwards  from 

a  value  of  At  to  a  value  for  |r,'2(f)|  to  a  value  for  |rg,(t)|.  Note  that  if  we  evaluate  (2.2-51)  at 
t  =  tm  +  T,then 

ro,i('m  +  t)|  =  ct - |r  [2{tm  +  r)|,  (2.2-52) 

and  upon  substituting  (2.2-49)  and  (2.2-50)  into  (2.2-52),  we  obtain 
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(2.2-53) 


ro..(^  +  T)|  =  |r;, 


where  jr^,  is  given  by  (2.2-32). 

Let  us  next  solve  for  the  unit  vectors  n'0l  and  n'u  so  that  we  can  evaluate  the  scattering 
function  g,(  /,nQ,,n,'2)  and,  hence,  the  complex  frequency  response  of  the  ocean.  Let  us  begin  with 
the  unit  vector  n'o  t .  From  (2.2-10), 


”o.i  =  r7"i(ro,i  + At'V,). 


(2.2-54) 


L  0,1 


Substituting  (2.2-30)  into  (2.2:54)  yields 


where 
we 


(2.2-55) 


r o,(r) |  is  given  by  (2.2-51),  r01  is  given  by  (2.1-6),  and  r  is  given  by  (2.2-49).  Note  that  if 
evaluate  (2.2-55)  at  t=tm  +  T,  then 


1 


ro,i  + 


ro,.(^  +  T) 


(2.2-56) 


|rw(tm  +  T)|[ 

and  upon  substituting  (2.2-53)  and  (2.2-30)  into  (2.2-56),  we  obtain 

<.('m +  *)  =  <„  (2.2-57) 

where  n'0]  is  given  by  (2.2-54).  And  by  referring  to  (2.1-35)  and  (2.1-36),  we  can  write  that 


3o,i(0  =  cos"  Wo/O,  t>tm  +  T, 


(2.2-58) 


and 


(2.2-59) 


where  u'0ft),  v'0l(t),  and  w^,(0  are  the  dimensionless,  time-varying  direction  cosines  with  respect 
to  the  X,  Y ,  and  Z  axes,  respectively,  associated  with  the  time- varying  unit  vector  n'0l(f)  given  by 
(2.2-55),  and  T  is  given  by  (2.2-49).  Equations  (2.2-58)  and  (2.2-59)  are  the  angles  of  incidence  at 
the  discrete  point  scatterer. 

From  (2.2-18), 
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1 


(2.2-60) 


«1,2=' 


1,2 


Substituting  (2.2-42)  into  (2.2-60)  yields 


(2.2-61) 


where  |r,'2(0|  is  given  by  (2.2-44),  r,  2  is  given  by  (2.1-14),  At  is  given  by  (2.2-48),  and  x  is  given 
by  (2.2-49).  Note  that  if  we  evaluate  (2.2-61)  at  t  =  tm  +  x,  then 


Ki(l + *)  = 


r'aitm  +  *) 


r.,2- 


T- 


+  ?)[ 


v, 


(2.2-62) 


and  upon  substituting  (2.2-49),  (2.2-50),  and  (2.2-30)  into  (2.2-62),  we  obtain 

*«('«  + =  (2-2-63) 

where  h[2  is  given  by  (2.2-60).  And  by  referring  to  (2.1-44)  and  (2.1-45),  we  can  write  that 


and 


0[2{t)  =  cos-1  w'i2(t) ,  t>  tm  +  r , 


(2.2-64) 


yr'2(r)  =  tan-1 


■\ 


J 


(2.2-65) 


where  v,'2(t),  and  w'l  2(t)  are  the  dimensionless,  time-varying  direction  cosines  with  respect 

to  the  X,  Y,  and  Z  axes,  respectively,  associated  with  the  time- varying  unit  vector  n[2{t)  given  by 

(2.2-61),  and  x  is  given  by  (2.2-49).  Equations  (2.2-64)  and  (2.2-65)  are  the  angles  of  scatter  at  the 
receiver. 

Note  that  since  the  unit  vectors  n'0l(t)  and  h[2{t)  given  by  (2.2-55)  and  (2.2-61),  respectively, 
are  functions  of  time,  the  scattering  function  is  also  a  function  of  time,  that  is, 

8lf,<vKi)=  t>tm  +  x.  (2.2-66) 

With  the  use  of  (2.2-22)  and  (2.2-23),  and  by  replacing  the  real  wavenumber  k  in  (2.2-23) 
with  the  complex  wavenumber  K  given  by  (2.1-20),  we  can  summarize  our  results  as  follows:  for 
the  bistatic  scattering  problem  shown  in  Fig.  2.2-1,  the  time-harmonic  velocity  potential  in  squared- 
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meters  per  second  incident  upon  the  receiver  at  r2  =  (x2,y2,z2),  due  to  a  time-harmonic  point  source 
at  r0  =  Cx0,y0,z0)  and  a  moving  discrete  point  scatterer  initially  at  r,  =  (jc,,y,,z,),  is  given  by 

yM (t,r2)  =  S0Hm  (t,r2  |/,r0)exp(+y2ff/r),  t>  tm  +  r,  (2.2-67) 


where  S0  is  the  source  strength  in  cubic  meters  per  second. 


(2.2-68) 


is  the  time-variant,  space-variant,  complex  frequency  response  of  the  ocean  at  frequency  /  hertz, 

s.(  /.<(*),««(■ ')) s  (2.2-69) 

is  the  scattering  function  of  the  discrete  point  scatterer  in  meters,  the  angles  of  incidence  0'Ql(t)  and 
Vo,i(t)  are  given  by  (2.2-58)  and  (2.2-59),  respectively,  the  angles  of  scatter  0f  2(O  and  \f/'U2(t)  are 
given  by  (2.2-64)  and  (2.2-65),  respectively,  cc(f)  is  the  real,  frequency-dependent,  attenuation 

coefficient  in  nepers  per  meter,  |i*Q,(r)|  is  given  by  (2.2-51),  |r'2(r)|  is  given  by  (2.2-44),  the  real 

wavenumber  k  in  radians  per  meter  is  given  by  (2.1-5),  and  the  time  delay  T  in  seconds  is  given  by 
(2.2-49). 

Let  us  next  relate  the  scattering  function,  the  differential  scattering  cross  section,  and  target 
strength  of  the  discrete  point  scatterer.  The  target  strength  (TS)  is  defined  as  follows  [6]: 


(2.2-71) 


is  the  differential  scattering  cross  section  with  units  of  squared  meters,  /  (r,'(0)  and  /  (r2)  are 

the  time-average,  incident  and  scattered  intensities,  respectively,  with  units  of  watts  per  squared 
meter,  is  the  scattering  function  of  the  discrete  point  scatterer  with  units  of 


26 


meters,  Aref  is  a  reference  cross-sectional  area  commonly  chosen  to  be  equal  to  1  m2 ,  and  the  time 
delay  r  in  seconds  is  given  by  (2.2-49). 


Example  2.2-1  Monostatic  Scattering  Geometry 

In  this  example  the  discrete  point  scatterer  is  in  motion  with  constant  velocity  vector  V,,  but 
the  scattering  geometry  is  monostatic  versus  bistatic.  For  a  monostatic  ( backscatter )  scattering 
geometry,  both  the  transmitter  and  receiver  are  located  at  the  same  position,  that  is, 


Substituting  (2.2-72)  into  (2.1-14)  yields 

r2  =  r0. 

(2.2-72) 

where  r0  I  is  given  by  (2.1-6).  Therefore, 

r..2=-ro,P 

(2.2-73) 

and 

rU2  =  r0A 

(2.2-74) 

A  A 

nh2  ~  ~n0,\  ' 

(2.2-75) 

And  upon  substituting  (2.1-73)  into  (2.2-18),  we  obtain 

ri,2  “  r0,l » 

(2.2-76) 

where  r^,  is  given  by  (2.2-10).  Therefore, 


and 


(2.2-77) 


"V  _  A/ 

ni,2  ~  ~  M0,P 


(2.2-78) 


and  upon  generalizing. 


<2(0|  =  |iw(0| 


(2.2-79) 


and 


<2  (*)  =  -<,(*)•  (2.2-80) 

With  the  use  of  (2.2-79)  and  (2.2-80),  the  time-variant,  space-variant,  complex  frequency 
response  of  the  ocean  given  by  (2.2-68)  reduces  to 
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hm (?>r2  |/>ro)  =  Si( /,»o,i(0-Wo,i(0)— ^  ^exp[-72/:|r;,(o|],  t>tm  +  t. 


(2.2-81) 


where 


-«o,,(0)  =  g,(  f,6'0J(t),y/'0l(t),ft  -  e'0l(t),7r  +  y/'0](t)). 

(2.2-82) 

0^,(0  and  y/'0l(t)  are  the  angles  of  incidence  given  by  (2.2-58)  and  (2.2-59),  respectively,  and 
0f,2(O  =  ft-  0q,i(O  and  W\, 2(0  =  ft  +  Woi(0  are  the  angles  of  scatter.  Substituting  (2.2-79)  into  (2.2- 
44),  (2.2-74)  and  (2.2-75)  into  (2.2-46)  and  (2.2-47),  and  (2.2-77)  into  (2.2-49)  yields 

|ro.i(0 

-S{t)  ±jg\t)  +  4s?(?(t) 

= - — - ,  t>tm  +  T, 

(2.2-83) 

where 

rf-.-rny. 

\  c  J 

(2.2-84) 

V 2  h  •  V 

g>(r)  =  2-^-Ar  +  2r0I-2^ - -, 

C  ’  c 

(2.2-85) 

m  =  Vj2(Ar)2  +  2  r0I(n01  •  V,)A  t  +  r2„ 

(2.2-86) 

A  t  =  t~tm,  t>tm  +  T, 

(2.2-87) 

and 


c 


(2.2-88) 


is  the  time  delay  in  seconds,  where  the  constant  value  of  range  ,  is  given  by  (2.2-32). 


2.3  All  Three  Platforms  In  Motion 

In  this  section  we  will  analyze  the  bistatic  scattering  problem  shown  in  Fig.  2.3-1.  All  three 
platforms  -  the  transmitter,  discrete  point  scatterer,  and  receiver  -  are  in  motion.  Motion  corresponds 
to  a  time-variant  problem.  As  mentioned  in  the  Introduction,  the  three  platforms  will  be  treated  as 
being  in  an  unbounded,  homogeneous  ocean  medium.  Although  the  propagation  of  sound  between 
the  source  and  discrete  point  scatterer,  and  between  the  discrete  point  scatterer  and  receiver  can  be 
treated  as  transmission  through  linear,  time-variant,  space-invariant  filters;  the  overall  solution  for 
this  bistatic  scattering  problem  corresponds  to  transmission  through  a  linear,  time-variant,  space- 
variant  filter.  The  presence  of  a  discrete  point  scatterer  in  an  unbounded,  homogeneous  fluid 
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medium  (i.e.,  a  fluid  medium  with  constant  speed  of  sound  and  ambient  density)  causes  the  medium 
to  be  space-variant. 

The  velocity  vectors  of  the  transmitter,  V0,  the  discrete  point  scatterer,  Vp  and  the  receiver,  V2, 
are  given  by 

V0  =  Vv  (2.3-1) 

V,  =  V^,  (2.3-2) 

and 

V2  =  V2«v  (2.3-3) 

where  V0,  V{,  and  V2  are  the  speeds  in  meters  per  second  of  the  transmitter,  the  discrete  point 
scatterer,  and  the  receiver,  respectively,  and  nVo ,  ,  and  are  the  dimensionless  unit  vectors  in 

the  directions  of  V0,  V,,  and  V2,  respectively.  The  velocity  vectors  given  by  (2.3-1)  through  (2.3-3) 
are  constant,  that  is,  the  speeds  and  directions  are  constants  -  there  is  no  acceleration.  Motion 
begins  at  time  t  =  tm  seconds.  We  will  model  the  propagation  of  sound  from  the  time  motion 
begins. 

Since  all  three  platforms  are  now  in  motion,  the  position  vectors  from  the  origin  to  the 
transmitter,  discrete  point  scatterer,  and  receiver  -  denoted  by  J£0(0,  /£,(0>  and  /R2(t),  respectively  - 
are  functions  of  time  given  by 


3>0(0  =  r0  +  A/V0, 

t>  t. 

m  ’ 

(2.3-4) 

£,(0  =  r,  +  ArV,, 

*■«*. 

IV 

3^ 

(2.3-5) 

3>2(r)  =  r2  +  ArV2, 

/>  rm, 

(2.3-6) 

where  r0  =  (x0,y0,z0),  r,  =  (*,, >,,<:,),  and  r2  =  (x2,y2,z2)  are  the  position  vectors  from  the  origin  to 
the  transmitter,  discrete  point  scatterer,  and  receiver,  respectively,  when  motion  begins  (see  Fig.  2.3- 
1),  and 


A  1=1-1.,  <><„■  (2.3-7) 

Note  that  R0(tm)  =  r0,  ^  ( rm )  =  r,,  and  %2(tm)  =  r2.  Instead  of  trying  to  solve  this  bistatic  scattering 
problem  directly  with  all  three  platforms  in  motion,  we  will  first  create  an  equivalent  problem 
involving  the  transmitter  and  the  discrete  point  scatterer  where  we  can  treat  the  transmitter  (sound 
source)  as  being  motionless.  We  will  then  create  a  second  equivalent  problem  involving  the  discrete 
point  scatterer  and  the  receiver  where  we  can  treat  the  discrete  point  scatterer  (acting  as  a  sound 
source)  as  being  motionless.  This  can  be  accomplished  by  working  with  relative  velocity  vectors. 

When  motion  begins  at  time  tm ,  the  scalar  component  of  V0  in  the  direction  of  V,  is  given  by 
(see  Fig.  2.3-2) 


Therefore,  the  velocity  vector  of  the  discrete  point  scatterer  relative  to  the  velocity  vector  of  the 
transmitter  in  the  direction  of  the  velocity  vector  of  the  discrete  point  scatterer  hv  is  given  by 


V,o  =  V1  -  («V,  *  VoK  =  [Vi-  Vo(nv,  •  «V0  )]«V 


(2.3-9) 


By  using  the  relative  velocity  vector  V,  0,  the  transmitter  (sound  source)  can  be  treated  as  being 
motionless.  Therefore,  when  motion  begins,  the  source  distribution  xM(t, r)  at  time  t  and  position 
r  =  (x,y,z)  will  be  treated  as  a  motionless,  time -harmonic,  point  source  with  units  of  inverse 
seconds,  that  is,  let 


xM(t,r)  =  S08(r-r0)exp(+j2nft),  (2.3-10) 

where  S0  is  the  source  strength  in  cubic  meters  per  second,  the  impulse  function  <5(r-r0),  with 
units  of  inverse  cubic  meters,  represents  a  point  source  at  r0  =  (*0,y0,z0),  a°d  /  is  frequency  in 
hertz.  The  sound  source  has  been  turned  on  forever,  that  is,  since  t  =  —°°.  In  addition,  since  we  are 
now  working  with  the  relative  velocity  vector  V,  0,  we  need  to  introduce  the  new  position  vector 

*r,(0  =  r1+AfVu,  t>tm,  (2.3-11) 

where  At  is  given  by  (2.3-7).  Compare  (2.3-11)  with  (2.3-5).  Note  that  £>1<1,0)(fm)  =  rl.  Also  note 
that  if  the  transmitter  is  not  in  motion,  then  Vo  =  0,  and  as  a  result,  V10  =  V,  [see  (2.3-9)]  and 

&,(,’0)(t)  =  (t)  [see  (2.3-1 1)  and  (2.3-5)].  And  if  the  discrete  point  scatterer  is  not  in  motion,  then 

V,  =  0  and  the  relative  velocity  vector  V,  0  given  by  (2.3-9)  is  undefined.  In  this  case  we  set 

V,  0  =  0,  and  as  a  result,  ^,<10)(0  =  %ft)  =  r, . 

When  the  transmitted  acoustic  field  is  first  incident  upon  the  discrete  point  scatterer  at  some 
time  t'  seconds  where  f  >  tm,  the  position  vector  from  the  origin  to  the  discrete  point  scatterer  is 
given  by  [see  (2.3-1 1)  and  Fig.  2.3-3] 

r,'=  *,(WV)  =  r,  +  At'  VI0,  f  >  tm,  (2.3-12) 


where 


^'=t'-tm,  t’>tm.  (2.3-13) 

The  propagation  of  sound  between  the  source  and  discrete  point  scatterer  can  be  modeled  as 
transmission  through  a  linear,  time-variant,  space -invariant  filter.  Therefore,  the  acoustic  field 
(velocity  potential)  incident  upon  the  discrete  point  scatterer  at  time  t'  and  position  r,' =  (x',  y'z.' ) 
is  given  by  [7] 


where 


y )  S0HM(f,T{  r0 |/jexp(+j27r/t  ),  t  >tm,  (2.3-14) 


i,\_  exp(-y*|r;-r0|) 

lJ) ““  47r|r'-r0| 


(2.3-15) 
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Figure  2.3-3.  Bistatic  scattering  geometry  when  the  transmitted  acoustic  field  is  first  incident 
upon  the  discrete  point  scatterer  at  time  t'  seconds  and  when  the  scattered  acoustic  field  is  first 
incident  upon  the  receiver  at  time  t  seconds  where  t>  t'  >  tm.  Point  0,  P0(r0),  is  the  transmitter; 

point  1,  /}(rj),  is  the  discrete  point  scatterer;  and  point  2,  P2(r2),  is  the  receiver.  All  three 
platforms  are  in  motion. 


is  the  time-variant,  space-invariant,  complex  frequency  response  of  the  ocean  at  frequency  /  hertz, 
and  k  is  the  wavenumber  in  radians  per  meter  given  by  (2.1-5)  and  is  repeated  below  for 
convenience: 


k  =  2  nf  jc  =  2  n/X . 


(2.1-5) 


By  referring  to  Fig.  2.3-3,  we  can  express  the  position  vector  from  the  point  source  to  the 
discrete  point  scatterer  at  time  t'  as 


ro,.  =  r.  -  r0’ 

and  upon  substituting  (2.3-12)  into  (2.3-16),  we  obtain 

t'  ~  r'  n'  —  r  4-  At'\ 
*0,1  'o.r'o.i  .  *o,iTiAt  vi,0’ 


(2.3-16) 


(2.3-17) 


where  r',= 


0,1 


and  is  repeated 


,  h'0l  is  the  dimensionless  unit  vector  in  the  direction  of  r'p  r0l  is  given  by  (2.1-6) 
below  for  convenience  (also  see  Fig.  2.3-3), 


ro,,  =  ri-ro> 

and  A t'  is  given  by  (2.3-13).  Therefore,  (2.3-14)  and  (2.3-15)  can  be  rewritten  as 


yM  (t',<)  =  S0Hm  (t',r'0,\f)exv(+j2nfO, 


t  >L, 


(2.1-6) 


(2.3-18) 


and 


exp) 

(-*!■ 

4  7t 

ro. 

l| 

(2.3-19) 


Let  us  now  create  a  similar  equivalent  problem  involving  the  discrete  point  scatterer  and  the 
receiver  where  we  can  treat  the  discrete  point  scatterer  (acting  as  a  sound  source)  as  being 
motionless.  When  the  transmitted  acoustic  field  is  first  incident  upon  the  discrete  point  scatterer  at 
time  t'  >  tm  seconds,  the  position  vector  from  the  origin  to  the  receiver  is  given  by  [see  (2.3-6)  and 
Fig.  2.3-3] 


r2'  =  *?2(0  =  r2  +  At'V2,  t'>tm,  (2.3-20) 

where  At'  is  given  by  (2.3-13).  Equation  (2.3-20)  indicates  that  after  At'  seconds,  the  receiver  - 
independent  of  the  transmitter  and  the  discrete  point  scatterer  -  travels  an  additional  distance  of 
A t'V2  meters  in  the  direction  nv  (see  Fig.  2.3-3).  Also  at  time  t',  the  scalar  component  of  V,  in 

the  direction  of  V2  is  given  by  (see  Fig.  2.3-4) 


•  V,  =  hv  •  V[nv  =  V{(hv  *nv). 


(2.3-21) 


Therefore,  the  velucily  vector  of  llie  receiver  relulive  to  the  velocity  vector  of  the  discrete  point 
scatterer  in  the  direction  of  the  velocity  vector  of  the  receiver  hVi  is  given  by 
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(2.3-22) 


By  using  the  relative  velocity  vector  V2 , ,  the  discrete  point  scatterer  (acting  as  a  sound  source)  can 
be  treated  as  being  motionless.  Therefore,  in  order  to  compute  the  acoustic  signal  incident  upon  the 
receiver,  we  treat  the  discrete  point  scatterer  at  time  t>t'>  tm  and  position  r,'  as  another 
motionless,  time-harmonic,  point  source  with  units  of  inverse  seconds,  that  is,  let  [see  (2.3-10)  and 
Fig.  2.3-3] 


x'M(t,  r)  =  S'<5(r  -  r;)exp(+j2nft),  •  (2.3-23) 

where  S'0  is  the  source  strength  in  cubic  meters  per  second,  and  the  impulse  function  <5(r-rj'), 
with  units  of  inverse  cubic  meters,  represents  a  point  source  at  r'  =  (x{,  y',  z\ ) ,  where  r,'  is  given  by 
(2.3-12).  The  source  strength  S0'  will  be  given  later. 

Since  we  are  now  working  with  the  relative  velocity  vector  V2 ,,  we  need  to  introduce  the  new 
position  vector 


^2’1) 


(0  =  r2  +  At"  V21, 


t>t'>t. 


m  ’ 


(2.3-24) 


where 


A  t"=f-t',  t>t'>tm,  (2.3-25) 

and  r2  is  given  by  (2.3-20).  Compare  (2.3-24)  with  (2.3-6).  Note  that  J t£2'x\t')  =  r2.  Also  note  that 
if  the  discrete  point  scatterer  is  not  in  motion,  then  V,  =  0,  and  as  a  result,  V2 ,  =  V2  [see  (2.3-22)] 
and  J522,l)(f)  =  &2(t)  for  t>t'>tm  [see  (2.3-24),  (2.3-20),  and  (2.3-6)].  And  if  the  receiver  is  not  in 
motion,  then  V2  =  0  and  the  relative  velocity  vector  V2 ,  given  by  (2.3-22)  is  undefined.  In  this  case 

we  set  V2 ,  =  0,  and  as  a  result,  %2'l\t)  =  %2(t)  =  r2  for  t>t'>tm. 

When  the  scattered  acoustic  field  is  first  incident  upon  the  receiver  at  some  time  t  seconds 
where  t>  t'  >  tm,  the  position  vector  from  the  origin  to  the  receiver  is  given  by  [see  (2.3-24)  and 
Fig.  2.3-3]  . 


r2"  =  ^2-1)(/)  =  r2'  +  A/"V2l,  t>t'>tm,  (2.3-26) 

where  A t"  is  given  by  (2.3-25).  Equation  (2.3-26)  indicates  that  after  A t"  seconds,  the  receiver 
travels  an  additional  distance  of  At"  V2  l  meters  in  the  direction  hv^  (see  Fig.  2.3-3).  By  referring  to 
Fig.  2.3-3,  it  can  be  seen  that 


rU  =  r2 


I  » 


(2.3-27) 


where  r"2  is  the  position  vector  between  the  discrete  point  scatterer  and  the  receiver  at  time 
t  >  t'  >  tm.  Substituting  (2.3-12)  and  (2.3-26)  into  (2.3-27)  yields 


ru=  rw  +  Al'(V2  - \0)  +  At"\2„ 


(2.3-28) 
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where  r"2  =  |r"2|,  n"2  is  the  dimensionless  unit  vector  in  the  direction  of  r"2,  r]2  is  given  by  (2.1- 
14)  and  is  repeated  below  for  convenience  (also  see  Fig.  2.3-3), 


ru-r2“rP  (2.M4) 

At'  is  given  by  (2.3-13),  and  At"  is  given  by  (2.3-25).  Equation  (2.3-28)  can  also  be  obtained 
directly  from  Fig.  2.3-3  -  without  the  use  of  (2.3-27)  -  simply  by  using  vector  addition.  Doing  so 
yields 


<2  =  V,,0  +  **1,2  +  A  t'\2  +  At" V2i, 

=  rIt2+A^(V2  -  V10)+Af"V21. 


(2.3-29) 


The  propagation  of  sound  between  the  discrete  point  scatterer  and  receiver  can  also  be 
modeled  as  transmission  through  a  linear,  time-variant,  space-invariant  filter.  Therefore,  the 
acoustic  field  (velocity  potential)  incident  upon  the  receiver  at  time  t  and  position  r2"  =  (x",y",z") 
due  to  a  point  source  at  r,'  =  (x'{,y'v  z[ )  is  given  by  [7] 

yM(t,r ")  =  S'0H  M(t,r" -r'\f)exp(+j2nft),  t>  t'  >  tm,  (2.3-30) 


where 


,,  I.  r»  r-|rt  expf-j*|r''- r,'  |) 


(2.3-31) 


With  the  use  of  (2.3-27),  (2.3-30)  and  (2.3-31)  can  be  rewritten  as 

yM & r2")  =  S'0Hm (r,r"|/)exp(+y2^/0,  t>t'>  tm, 

and 


(2.3-32) 


(2.3-33) 


The  source  strength  Sq  is  given  by  [see  (2.3-18)] 


S'o  =  So»»('X,k)*.(  /•<,.%).  P-3’34) 

where  /,n',,n"2)  is  the  scattering  function  of  the  discrete  point  scatterer  with  units  of  meters  as 
discussed  in  Section  2.1.  Since  the  transmitted  acoustic  field  is  first  incident  upon  the  discrete  point 
scatterer  at  time  t' ,  the  direction  of  wave  propagation  from  the  source  to  the  scatterer  is  given  by  the 
dimensionless  unit  vector  n0A  (see  Fig.  2.3-3).  Similarly,  since  the  discrete  point  scatterer  is  being 
treated  as  another  point  source  at  time  t'  and  position  r' ,  and  the  scattered  acoustic  field  is  first 
incident  upon  the  receiver  at  time  t  >  t'  >  tm,  the  direction  of  wave  propagation  from  the  scatterer  to 
the  receiver  is  given  by  the  dimensionless  unit  vector  n"2  (see  Fig.  2.3-3).  Later  in  this  section,  we 
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will  show  that  the  scattering  function  is  also  a  function  of  time  because  the  unit  vectors  are  actually 
time-varying.  We  will  then  show  how  to  express  the  scattering  function  as  a  function  of  frequency 
and  time-varying  angles  of  incidence  and  scatter  instead  of  time-varying  unit  vectors.  The  use  of 
unit  vectors  is  meant  as  a  shorthand  notation. 

Let  us  now  begin  the  process  of  obtaining  a  final  expression  for  the  time-harmonic  velocity 
potential  incident  upon  the  receiver.  Substituting  (2.3-34)  into  (2.3-32)  yields 

(' 

yM  <Lr")  =  S0Hm  {t\r'a\f)g(f,K.vK2)HM  (t,r"2\f)e\p(+j2nft),  t>t'>tm, 

(2.3-35) 

or,  equivalently, 

yM  ri )  =  So tiM  [t, r2"  | f,r0)e\p{+j2nft),  t>t'>tm,  (2.3-36) 


where 


H\i  (*» r"  I  /» ro ) =  Hm  Fq,  |  /)  g ,  ( /,  «o,p  n",2  r"2|  /  j 

=  8i(f*K.vKi) 

=  8i{f’KiX,2) 


exp 

~jk( 

ro,i 

+ 

c)' 

167T2 

r0.t 

ru 

exp 

r,-r0  +  Af/Vli0 

+ 

r2-r,+A('(Vi-Vl,0)+A<"V1,|)] 

16tt2 

Tj  —  r0  +  A t  Vli0 

r2 

-  r,  +  At\  V2  -  V,  0 ) + At"  V2 , 

(2.3-37) 


At'  is  given  by  (2.3-13),  V,  0  is  given  by  (2.3-9),  At"  is  given  by  (2.3-25),  V2 ,  is  given  by  (2.3- 
22),  and 


or 


t  =  t'  + 


(2.3-38) 


(2.3-39) 


where  r,'2  is  given  by  (2.3-28).  Note  that  if  the  bistatic  scattering  problem  shown  in  Fig.  2.3-1 
corresponded  to  transmission  through  a  space-invariant  filter,  then  the  complex  frequency  response 
given  by  (2.3-37)  would  be  a  function  of  the  vector  spatial  difference  r"  -  r0 ,  which  it  is  not  [  r"  is 
given  by  (2.3-26)]. 
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In  order  to  evaluate  the  complex  frequency  response  given  by  (2.3-37),  we  must  derive 


solutions  for  the  unit  vectors  h'0i  and  n"2,  and  the  ranges 
quantities.  Let  us  begin  with  |r„  ,1 .  Since 


■  0,1 


and 


■1,2 


in  terms  of  known 


ro,.  =  ro,.#ro,i> 


(2.3-40) 


substituting  (2.3-17)  into  (2.3-40)  yields 


<,  =  (r0 1  +  At '  V,  0 )  •  (r0 ,  +  At'  V,  0 ) . 


(2.3-41) 


Expanding  the  right-hand  side  of  (2.3-41)  and  taking  the  square  root  of  both  sides  of  the  resulting 
equation  yields 


r0,l  —  r0,l 


1  +  —  («o.,*VI,0)At'  + 


'0,1 


rYuAA2 


1/2 


(2.3-42) 


or 


r'  1  —  r  . 

ro,i  |  —  V 

v? 

N) 

h 

H _ 

(«0,1*«V1)  + 


1  VWA  f 


2  r, 


0,1 


1/2 


(2.3-43) 


Although  (2.3-42)  and  (2.3-43)  are  exact  expressions  for  the  range  |r',|,  we  cannot  compute 

a  value  for  |r„  ,|  until  we  derive  an  expression  for  A t'  in  terms  of  |r^  ,| .  By  referring  to  Fig.  2.3-3,  it 
can  be  seen  that 


> 

II 

O  n, 

t'>t. 

m 

c 

(2.3-44) 


Substituting  (2.3-44)  into  (2.3-42)  and  squaring  both  sides  of  the  resulting  equation  yields  the 
following  second-order  polynomial 


A 


(2.3-45) 


with  exact  solution 


Ir' 

B±^B2  +  4AC 

|ro,i 

2  A 

(2.3-46) 
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where 


and 


B  =  2r0l  —  -  ^',0 
c 


(2.3-47) 

(2.3-48) 


C=r2or  (2.3-49) 

The  solution  given  by  (2.3-46)  is  the  constant  value  of  range  between  the  point  source  and  the 
discrete  point  scatterer  when  the  transmitted  acoustic  field  is  first  incident  upon  the  discrete  point 
scatterer  at  time  instant  t'  after  motion  begins  at  time  instant  tm  where  f  >  tm.  The  decision  to  use 
either  the  plus  or  minus  sign  in  (2.3-46)  is  dictated  by  the  fact  that  range  must  be  positive.  Let  us 
solve  for  r"2  next. 

Since 


substituting  (2.3-29)  into  (2.3-50)  yields 

|ru|2  =  [ri,2  +  Af'( V2  -  V, 0) + At" V2 ,] •  [r,  2  + Ar'(  V2  -  V, 0)  +  Ar" V2 ,].  (2.3-51) 


Expanding  the  right-hand  side  of  (2.3-51)  and  taking  the  square  root  of  both  sides  of  the  resulting 
equation  yields 


(2.3-52) 


or 
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(2.3-53) 

Although  (2.3-52)  and  (2.3-53)  are  exact  expressions  for  the  range  |r"2|,  we  cannot  compute 
a  value  for  |r"2|  until  we  derive  appropriate  expressions  for  A t' ,  A t",  and  At'At" .  We  will  first 
solve  for  a  constant  value  for  |r"2|  by  expressing  At' ,  At",  and  At'At"  in  terms  of  the  known 
constant  range  ,|  [see  (2.3-46)]  and  |r"2|.  We  will  then  solve  for  a  time-varying  |r"2|  by 
expressing  At',  At",  and  At'At"  exclusively  in  terms  of  |r"2|.  Since  (2.3-44)  already  expresses 

At'  in  terms  of  ,1,  we  begin  by  substituting  (2.3-38)  into  (2.3-25)  which  yields  (also  see  Fig. 
2.3-3) 

(2.3-54) 

and  upon  multiplying  (2.3-44)  by  (2.3-54),  we  obtain 

(2.3-55) 

The  solution  for  |r"2|  can  be  obtained  by  substituting  (2.3-54)  into  (2.3-52)  and  squaring 
both  sides  of  the  resulting  equation.  Doing  so  yields  the  following  second-order  polynomial 

=  °’  (2.3-56) 

with  exact  solution 

(2.3  57) 
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where 


4  -  l _ i 

^0  ~ 1  2  1 

C 

(2.3-58) 

A,  +  A.At' 

*•-  c  ’ 

(2.3-59) 

£0=A2(Ar')2  +  A,A/'  +  4 

(2.3-60) 

^1  =  ^  *1,2  [^1,2  •  (V2  —  V|  0)], 

(2.3-61) 

^=|v,-vj2. 

(2.3-62) 

A3  =  2r12  («1|2*  V2,i), 

(2.3-63) 

ii 

(2.3-64) 

Aj=2[Vu.(V2-Vu)], 

(2.3-65) 

A  t'=t'-t  = 


1 0,1 


f>tm. 


(2.3-66) 


where  ,|  is  given  by  (2.3-46).  The  solution  given  by  (2.3-57)  is  the  constant  value  of  range 

between  the  discrete  point  scatterer  and  the  receiver  when  the  scattered  acoustic  field  is  first  incident 
upon  the  receiver  at  time  instant  t  after  motion  begins  at  time  instant  tm  where  t>  t'  >  tm.  The 
decision  to  use  either  the  plus  or  minus  sign  in  (2.3-57)  is  dictated  by  the  fact  that  range  must  be 
positive. 

As  we  previously  mentioned,  we  will  now  solve  for  a  time-varying  |r"2|  by  expressing  At', 

At",  and  At' At"  exclusively  in  terms  of  |r"2|.  We  begin  by  substituting  (2.3-39)  into  (2.3-13) 
which  yields 


A  t'=t-tm 

m 


t>l'>t 


and  upon  substituing  (2.3-7)  into  (2.3-67),  we  obtain 


At' =  A  t- 


t>t'>t. 

m 


(2.3-67) 


(2.3-68) 


And  upon  multiplying  (2.3-68)  by  (2.3-54), 
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(2.3-69) 

The  solution  for  |r"2|  can  be  obtained  by  substituting  (2.3-68),  (2.3-54),  and  (2.3-69)  into 

(2.3-52)  and  squaring  both  sides  of  the  resulting  equation.  Doing  so  yields  the  following  second- 
order  polynomial 

^|r"  (t)f  +  (f)|r£(f)|  -  m  =  0,  t>tm  +  t,  (2.3-70) 

with  exact  solution 

(2.3-71) 


where 

(A2  +  A4  —  A5) 

^  =  1 - 2 - . 

c 

(2A2-A5)Ar  +  (A.-A3) 

^(0  =  - — 1 — - — — -, 

c 

<?(0  =  A2(A02+A,Ar  +  r22, 

A{  through  A5  are  given  by  (2.3-61)  through  (2.3-65), 

A  t>tm  +  r, 

and  (see  Fig.  2.3-3) 


(2.3-76) 

is  the  time  delay  in  seconds  (the  amount  of  time  it  takes  for  the  transmitted  acoustic  signal  to  begin 
to  appear  at  the  receiver  after  motion  begins  at  time  instant  tm )  where  the  constant  values  of  range 

| Tq  ! |  and  r"2|  are  given  by  (2.3-46)  and  (2.3-57),  respectively.  It  is  important  to  note  that  if  (2.3- 

70)  is  evaluated  at  t  =  tm  +  t,  then  it  can  be  shown  that  (2.3-70)  can  be  rewritten  in  the  form  of  (2.3- 
56)  and,  as  a  result, 

ru«„+T)|  =  |>4  (2-3-77) 


(2.3-72) 

(2.3-73) 

(2.3-74) 

(2.3-75) 
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where  |r"2|  is  given  by  (2.3-57). 

Now  that  we  have  an  exact  solution  for  |r"2(o|  as  given  by  (2.3-71),  we  can  use  it  to  obtain  an 

exact  solution  for  the  time-varying  range  |ro,(r)|  as  follows.  With  the  use  of  (2.3-44)  and  (2.3-68), 
we  can  write  that 

(2.3-78) 

where  A t  is  given  by  (2.3-75),  |r"2(/)|  is  given  by  (2.3-71),  and  r  is  given  by  (2.3-76).  Equation 
(2.3-78)  indicates  that  if  we  are  given  a  value  of  A t,  which  determines  the  value  of  |r"2(i)|,  then  we 
can  use  those  two  values  to  compute  what  |ro  ,(r)|  must  have  been  -  we  are  working  backwards  from 

a  value  of  At  to  a  value  for  |r"2(r)|  to  a  value  for  |ioj(f)|.  Note  that  if  we  evaluate  (2.3-78)  at 
t  =  t„  +  t,  then 

tit  7 

|ro,.('m  +  T)|  =  cT-|r"  (tm  +  t)|,  (2.3-79) 

and  upon  substituting  (2.3-76)  and  (2.3-77)  into  (2.3-79),  we  obtain 

ro.i(*m  +  T)|  =  |ro,i|’  (2.3-80) 

where  |r„  ,|  is  given  by  (2.3-46). 

Let  us  next  solve  for  the  unit  vectors  n'0l  and  n"2  so  that  we  can  evaluate  the  scattering 
function  /, n'0A,n"2)  and,  hence,  the  complex  frequency  response  of  the  ocean.  Let  us  begin  with 
the  unit  vector  n'ox .  From  (2.3-17), 

«;.1=i4|(r»j+A,'v'.«)-  <2-3-81) 

|r0,l| 

Substituting  (2.3-66)  into  (2.3-81)  yields 

(2.3-82) 

where  1^,(0]  is  given  by  (2.3-78),  r01  is  given  by  (2.1-6),  V10  is  given  by  (2.3-9),  and  T  is  given 
by  (2.3-76).  Note  that  if  we  evaluate  (2.3-82)  at  t  =  tm  +  t,  then 


«o,.('m  +  T)  = 


1 

r0,.('m  +  T)|,r 

r0,l(^m  +  T) 

r0,l  + 

M,0 

c  J 

(2.3-83) 


and  upon  substituting  (2.3-80)  and  (2.3-44)  into  (2.3-83),  we  obtain 

(23-84> 

where  n'0l  is  given  by  (2.3-81).  And  by  referring  to  (2.2-58)  and  (2.2-59),  we  can  write  that 


0o  ft)  =  cos-1  w'(t),  t>tm  +  T, 


'o,iv 


(2.3-85) 


and 


(2.3-86) 


where  u'Ql(t),  v'0ft),  and  w'0ft)  are  the  dimensionless,  time-varying  direction  cosines  with  respect 
to  the  X,  Y,  and  Z  axes,  respectively,  associated  with  the  time-varying  unit  vector  n'Ql(t)  given  by 
(2.3-82),  and  r  is  given  by  (2.3-76).  Equations  (2.3-85)  and  (2.3-86)  are  the  angles  of  incidence  at 
the  discrete  point  scatterer. 

From  (2.3-28), 


«s= rril^ + A,'(  Vj  ■  v'-°) + AI"  v4 


(2.3-87) 


1 1,2 1 


Substituting  (2.3-68)  and  (2.3-54)  into  (2.3-87)  yields 


(2.3-88) 


r"2(o|  is  given  by  (2.3-71),  rl2  is  given  by  (2.1-14),  A t  is  given  by  (2.3-75),  V10  is  given 
-9),  V21  is  given  by  (2.3-22),  and  t  is  given  by  (2.3-76).  Note  that  if  we  evaluate  (2.3-88)  at 


where 
by  (2.3 
t  =  tm  +  T,  then 


n"jL  +  t)  = 


ri,2  + 


T- 


4  T)| 


r  "(tm  +  t)| 

(V2- vli0)  +  J-;-— ^ — [V2,I 


(2.3-89) 
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and  upon  substituting  (2.3-76),  (2.3-77),  (2.3-44),  and  (2.3-54)  into  (2.3-89),  we  obtain 

nUtm  +  r)  =  h['2,  (2.3-90) 

where  h"2  is  given  by  (2.3-87).  And  by  referring  to  (2.2-64)  and  (2.2-65),  we  can  write  that 

(2.3-91) 
and 

(2.3-92) 

where  u"2(t),  v"2(t),  and  w"2(t)  are  the  dimensionless,  time-varying  direction  cosines  with  respect 
to  the  X ,  Y ,  and  Z  axes,  respectively,  associated  with  the  time- varying  unit  vector  h['2{t)  given  by 
(2.3-88),  and  x  is  given  by  (2.3-76).  Equations  (2.3-91)  and  (2.3-92)  are  the  angles  of  scatter  at  the 
receiver. 

Note  that  since  the  unit  vectors  n'0l(t)  and  h['2{i)  given  by  (2.3-82)  and  (2.3-88),  respectively, 
are  functions  of  time,  the  scattering  function  is  also  a  function  of  time,  that  is, 

•  t>tm  +  r.  (2.3-93) 

Let  us  next  derive  an  expression  for  the  time- varying  position  vector  to  the  receiver  when  the 
scattered  acoustic  field  is  incident  upon  the  receiver.  We  begin  by  substituting  (2.3-20)  into  (2.3- 
26).  Doing  so  yields 

r"  =  r2  +  A?'  V2  +  At"  V2I ,  (2.3-94) 

which  is  the  position  vector  from  the  origin  to  the  receiver  when  the  scattered  acoustic  field  is  first 
incident  upon  the  receiver  at  time  instant  t  after  motion  begins  at  time  instant  tm  where  t>  t'  >  tm. 
And  upon  substituting  (2.3-68)  and  (2.3-54)  into  (2.3-94),  we  obtain 


(2.3-95) 

where  A t  is  given  by  (2.3-75),  r"2(o|  is  given  by  (2.3-71),  V2,  is  given  by  (2.3-22),  and  x  is  given 
by  (2.3-76).  Note  that  if  we  evaluate  (2.3-95)  at  t=tm  +  x,  then 


and  upon  substituting  (2.3-76)  and  (2.3-77)  into  (2.3-96),  we  obtain 


r2'U,  +  t)  =  r2  + 


(2.3-97) 


where  ,|  and  | r"2 1  are  given  by  (2.3-46)  and  (2.3-57),  respectively.  If  we  further  substitute  (2.3- 
44)  and  (2.3-54)  into  (2.3-97),  then 


r?(t.  +  0  =  r2",  (2.3-98) 

where  r"  is  given  by  (2.3-94). 

With  the  use  of  (2.3-36)  and  (2.3-37),.  and  by  replacing  the  real  wavenumber  k  in  (2.3-37) 
with  the  complex  wavenumber  K  given  by  (2.1-20),  we  can  summarize  our  results  as  follows:  for 
the  bistatic  scattering  problem  shown  in  Fig.  2.3-1,  the  time-harmonic  velocity  potential  in  squared- 
meters  per  second  incident  upon  the  receiver  at  r"(f)  =  (x"(t),y"(t),z"(t)),  due  to  a  moving  time- 
harmonic  point  source  initially  at  r0  =  (x0,y0,z0)  and  a  moving  discrete  point  scatterer  initially  at 
r,  =(xvyl,zl),  is  given  by 


yM  (*>r2"(0)  =  S0Hm  (t,r"(t)\f,r0}exp(+j2nft),  t  >  tm  + T, 


(2.3-99) 


where  r  "(t)  is  given  by  (2.3-95),  50  is  the  source  strength  in  cubic  meters  per  second, 


exp 

rr  -f  »  —  rr  (  f  /f\ 

+ll 

£«>|) 

2.oV*J_  aYL'  r/J -i_L"/'rtlYI 

w| —  £i\  /  »**o,lw»^l,2w/ 

16tt2|i 

•wol 

(2.3-100) 


is  the  time- variant,  space- variant,  complex  frequency  response  of  the  ocean  at  frequency  /  hertz, 


(2.3-101) 


is  the  scattering  function  of  the  discrete  point  scatterer  in  meters,  the  angles  of  incidence  6'0  i(t)  and 
\ff'Q](t)  are  given  by  (2.3-85)  and  (2.3-86),  respectively,  the  angles  of  scatter  6"2{t)  and  V'uCO  316 
given  by  (2.3-91)  and  (2.3-92),  respectively,  a(/)  is  the  real,  frequency-dependent,  attenuation 

coefficient  in  nepers  per  meter,  jr^,(o|  is  given  by  (2.3-78),  |t"2(t)|  is  given  by  (2.3-71),  the  real 

wavenumber  k  in  radians  per  meter  is  given  by  (2.1-5),  and  the  time  delay  t  in  seconds  is  given  by 
(2.3-76). 

Let  us  next  relate  the  scattering  function,  the  differential  scattering  cross  section,  and  target 
strength  of  the.  discrete  point  scatterer.  The  target  strength  (TS)  is  defined  as  follows  [6]: 
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(2.3-103) 


is  the  differential  scattering  cross  section  with  units  of  squared  meters,  Iavg  (r,'(t))  and  /  (r2"(if)) 

are  the  time-average,  incident  and  scattered  intensities,  respectively,  with  units  of  watts  per  squared 
meter,  g,(  f  ,n'0l(t),h"2{t)}  is  the  scattering  function  of  the  discrete  point  scatterer  with  units  of 

meters,  Att(  is  a  reference  cross-sectional  area  commonly  chosen  to  be  equal  to  1  m2 ,  and  the  time 
delay  t  in  seconds  is  given  by  (2.3-76). 


Example  2.3-1  Simulation  of  Sections  2.1  and  2.2 

In  this  example  we  will  show  that  the  exact  results  derived  in  this  section  will  reduce  to  the 
exact  results  derived  in  Sections  2.1  and  2.2  when  appropriate  values  are  used  for  the  various 
parameters.  Let  us  begin  with  the  problem  discussed  in  Section  2.2. 

Simulation  of  Section  2.2 
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A  =  1- 


B  =  2r, 


Y 

c 

yv, 


0.1 


and 


C=r, 


0,1  • 


(2.3-108) 

(2.3-109) 

(2.3-110) 


Equations  (2.3-108)  through  (2.3-1 10)  are  identical  with  (2.2-33)  through  (2.2-35).  Therefore,  with 
the  use  of  (2.3-108)  through  (2.3-1 10),  the  exact  solution  for  the  constant  value  of  range  IrjJ  given 


by  (2.3-46)  reduces  to  the  exact  solution  given  by  (2.2-32).  Recall  that  r^,  is  the  constant  value  of 

range  between  the  point  source  and  the  discrete  point  scatterer  when  the  transmitted  acoustic  field  is 
first  incident  upon  the  discrete  point  scatterer  at  time  instant  t'  after  motion  begins  at  time  instant 
tm  where  t'  >  tm. 

m  m 

Substituting  V2=  0  and  (2.3-104)  through  (2.3-107)  into  (2.3-61)  through  (2.3-65)  yields 


and 


At  =  —  2r12(w12*  Vx>, 

(2.3-111) 

a2=v,2, 

(2.3-112) 

a3= 0, 

(2.3-113) 

•ti¬ 

ll 

© 

(2.3-114) 

© 

II 

(2.3-115) 

Substituting  (2.3-113)  through  (2.3-115)  into  (2.3-58)  and  (2.3-59),  and  repeating  (2.3-60)  for 
convenience  yields 


/*0=1,  (2.3-116) 

^0=0,  (2.3-117) 

and 

£0=A2(A/')2  +  A,A*'  +  /£l.  (2.3-118) 

Substituting  (2.3-1 16)  through  (2.3-1 18)  into  (2.3-57)  yields 

|r£|=^A2(A/')2  +  AlAr'  +  £,  (2.3-119) 

and  upon  substituting  (2.3-66),  (2.3-1 1 1),  and  (2.3-1 12)  into  (2.3-1 19),  we  finally  obtain 
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(2.3-120) 


where  |r,'2|  is  given  by  (2.2-40).  Therefore,  the  exact  solution  for  the  constant  value  of  range  |r"2| 
given  by  (2.3-57)  reduces  to  the  exact  solution  for  the  constant  value  of  range  |r,'2|  given  by  (2.2- 

40).  Recall  that  both  r"2  and  j  r 2 1  are  the  constant  values  of  range  between  the  discrete  point 

scatterer  and  the  receiver  when  the  scattered  acoustic  field  is  first  incident  upon  the  receiver  at  time 
instant  t  after  motion  begins  at  time  instant  tm  where  t>  t'  >tm. 

Substituting  (2.3-1 1 1)  through  (2.3-1 15)  into  (2.3-72)  through  (2.3-74)  yields 


(2.3-121) 


and 


V2 

S(t) -  2— A/-2r12 
c 


*L2»V, 

C 


(2.3-122) 


m  =  Vl2m2-2rl2(nl2»\l)M+r22,  (2.3-123) 

where  At  is  given  by  (2.3-75).  And  because  of  (2.3-120),  At  given  by  (2.3-75)  is  identical  to  At 
given  by  (2.2-48)  since  the  time  delay  r  given  by  (2.3-76)  is  identical  to  the  time  delay  r  given  by 
(2.2-49).  As  a  result,  (2.3-121)  through  (2.3-123)  are  identical  with  (2.2-45)  through  (2.2-47). 
Therefore,  with  the  use  of  (2.3-121)  through  (2.3-123),  the  exact  solution  for  the  time-varying  range 

|r"2(t)  given  by  (2.3-71)  reduces  to  the  exact  solution  for  the  time-varying  range  |r,'2(t)|  given  by 
(2.2-44),  that  is, 

|r"(()|  =  |r'2(0|.  (2.3-124) 

With  the  use  of  (2.3-124),  the  exact  solution  for  the  time-varying  range  |r^|(r)l  given  by  (2.3- 
78)  reduces  to  the  exact  solution  given  by  (2.2-51).  With  the  use  of  (2.3-104),  the  exact  solution  for 
the  time-varying  unit  vector  n'0A(t)  given  by  (2.3-82)  reduces  to  the  exact  solution  given  by  (2.2- 
55).  With  the  use  of  V2  =  0,  (2.3-104),  (2.3-106),  and  (2.3-124),  the  exact  solution  for  the  time- 
varying  unit  vector  n"2(t)  given  by  (2.3-88)  reduces  to  the  exact  solution  for  the  time-varying  unit 
vector  n[2{t)  given  by  (2.2-61),  that  is, 

n''2(t)  =  n[2(t).  (2.3-125) 

And  finally,  substituting  V2=  0  and  (2.3-106)  into  (2.3-95)  yields 

r"(r)  =  r2.  (2.3-126) 

Let  us  next  consider  the  problem  discussed  in  Section  2. 1 
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Simulation  of  Section  2.1 

In  Section  2.1  none  of  the  platforms  were  in  motion.  In  order  to  simulate  this  problem,  all  we 
need  to  do  for  the  most  part  is  to  set  V,  =  0  and  Vj  =  0  into  the  equations  already  derived  in  this 
example. 

Substituting  V,=  0  into  (2.3-108)  and  V,  =  0  into  (2.3-109),  and  repeating  (2.3-110)  for 
convenience  yields 


A  =  l, 

(2.3-127) 

B  =  0, 

.  (2.3-128) 

and 

o 

II 

pY> 

(2.3-129) 

Substituting  (2.3-127)  through  (2.3-129)  into  (2.3-46)  yields 

r0.l|=  r0.1  =  r0,l 

(2.3-130) 

as  the  constant  value  of  range  between  the  point  source  and  the  discrete  point  scatterer  [see  (2.1- 
26)]. 

Substituting  V,  =  0  into  (2.3-111)  and  V,  =  0  into  (2.3-112),  and  repeating  (2.3-113)  through 
(2.3- 115)  for  convenience  yields 

A,=  0, 

(2.3-131) 

A2=  0, 

(2.3-132) 

A3  =  0, 

(2.3-133) 

4*. 

II 

O 

(2.3-134) 

and 

d 

II 

(2.3-135) 

Substituting  (2.3-131)  and  (2.3-132)  into  (2.3-118)  and  (2.3-119),  and  repeating  (2.3-116)  and 
(2.3- 117)  for  convenience  yields 

* 

(2.3-136) 

S0=0, 

(2.3-137) 

ri,2’ 

(2.3-138) 

and 
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(2.3-139) 


=  r, 


1,2 


as  the  constant  value  of  range  between  the  discrete  point  scatterer  and  the  receiver  [see  (2.1-27)]. 
Note  that  with  the  use  of  (2.3-130)  and  (2.3-139),  the  time  delay  r  given  by  (2.3-76)  reduces  to  the 
time  delay  r  given  by  (2.1-25). 

Substituting  V,  =  0  and  Vl  =  0  into  (2.3-121)  through  (2.3-123)  yields 


*= 1, 

(2.3-140) 

m=o, 

(2.3-141) 

and 

m 

(2.3-142) 

Substituting  (2.3-140)  through  (2.3-142)  into  (2.3-71)  yields 

|ri,2(*)|  -  |ri,2  —  ri,2- 

(2.3-143) 

Equation  (2.3-143)  indicates  that  the  range  between  the  discrete  point  scatterer  and  the  receiver  does 
not  change  value  as  a  function  of  time  -  it  is  equal  to  the  constant  value  rl2  -  which  makes  physical 
sense  since  none  of  the  platforms  are  in  motion. 

When  none  of  the  platforms  are  in  motion,  the  time  delay  between  the  transmitter,  discrete 
point  scatterer,  and  receiver  is  constant,  that  is,  At  =  r,  where  T  is  given  by  (2.3-76).  Therefore,  with 
the  use  of  (2.3-76),  (2.3-130),  and  (2.3-139), 


A  t  =  T  = 


c  c 


Substituting  (2.3-143)  and  (2.3-144)  into  (2.3-78)  yields 


(2.3-144) 


(2.3-145) 


Equation  (2.3-145)  indicates  that  the  range  between  the  transmitter  and  the  discrete  point  scatterer 
does  not  change  value  as  a  function  of  time  -  it  is  equal  to  the  constant  value  r01  -  which  makes 
physical  sense  since  none  of  the  platforms  are  in  motion. 

Substituting  V,=  0  into  (2.3-104)  yields 


V,.o=0.  (2.3-146) 

With  the  use  of  (2.3-145)  and  (2.3-146),  the  exact  solution  for  the  time-varying  unit  vector  n'0{(t) 
given  by  (2.3-82)  reduces  to  the  exact  solution  for  the  constant  unit  vector  n01  given  by  (2.1-18), 
that  is, 


«O,l(0  =«(,..= 


0,1 


(2.3-147) 
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And  finally,  with  the  use  of  V2=  0,  (2.3-106),  (2.3-143),  and  (2.3-146),  the  exact  solution  for  the 
time-varying  unit  vector  n"2(t)  given  by  (2.3-88)  reduces  to  the  exact  solution  for  the  constant  unit 
vector  n12  given  by  (2.1-19),  that  is, 


n"(0=n,2  = 


1,2 


1.2 


(2.3-148) 


Example  2.3-2  Monostatic  Scattering  Geometry 

In  this  example  all  three  platforms  are  in  motion  with  constant  velocity  vectors  V0,  V,,  and 
V2,  but  the  scattering  geometry  is  monostatic  versus  bistatic.  For  a  monostatic  ( backscatter ) 
scattering  geometry,  both  the  transmitter  and  receiver  are  located  at  the  same  position,  that  is, 

r2  =  r0.  (2.3-149) 

Substituting  (2.3-149)  into  (2.1-14)  yields 

ri,2  =  -ro,P  (2.3-150) 


where  r0  ,  is  given  by  (2.1-6).  Therefore, 


and 


ria  ~  ro,i 


(2.3-151) 


«1,2  =  -«0.|- 

(2.3-152) 

And  since  both  the  transmitter  and  receiver  are  on  the  same  platform, 

II 

o< 

(2.3-153) 

II 

(2.3-154) 

and 

fly  — *  fly  • 
v2 

(2.3-155) 

Substituting  (2.3-153)  through  (2.3-155)  into  (2.3-22)  yields  the  following  expression  for  the 
velocity  vector  of  the  receiver  relative  to  the  velocity  vector  of  the  discrete  point  scatterer  in  the 
direction  of  the  velocity  vector  of  the  receiver  hVi  =  hVo : 

v2,  =  V0  -  («v, .  V,  )nVi  =  [  V„-  V,  (nv> .  \  )]nVj .  (2.3-156) 

The  time  variant,  space  variant,  complex  frequency  rc3pon3c  of  the  ocean  for  a  mono3tatic 

scattering  geometry  with  all  three  platforms  in  motion  is  given  by  (2.3-100)  in  conjunction  with 
(2.3-149)  through  (2.3-156). 
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Example  2.3-3  Synthetic  Aperture  Sonar  (SAS) 

In  this  example  we  will  consider  a  bistatic  scattering  geometry  consistent  with  a  synthetic 
aperture  sonar  (SAS)  trying  to  image  a  nonmoving  target  on  the  ocean  bottom.  Although  the 
transmitter  and  receiver  are  on  the  same  platform,  they  are  not  located  at  the  same  position,  that  is, 


r2*r0. 

(2.3-157) 

However,  (2.3-153)  through  (2.3-156)  are  applicable  in  this  example. 

As  was  mentioned,  the  discrete  point  scatterer  (target)  is  not  in  motion.  Therefore, 

v,=o. 

(2.3-158) 

and,  as  a  result,  the  relative  velocity  vector  V,  0  given  by  (2.3-9)  is  undefined.  Therefore, 

v,0=o 

(2.3-159) 

1 

Q- 

© 

II 

o _ 

II 

o 

(2.3-160) 

And,  upon  substituting  (2.3-158)  into  (2.3-156),  we  obtain 

V  =  V 

T  2.1  T  0 

(2.3-161) 

and 

v^KHv.l 

(2.3-162) 

-Vo- 

Let  us  solve  for  the  constant  value  of  range  |r(J  next.  Substituting  (2.3-160)  into  (2.3-47) 

and  (2.3-159)  into  (2.3-48),  and  repeating  (2.3-49)  for  convenience  yields 

A  =  1, 

(2.3-163) 

5  =  0, 

(2.3-164) 

and 

• 

C=rl. 

(2.3-165) 

Substituting  (2.3-163)  through  (2.3-165)  into  (2.3-46)  yields 


(2.3-166) 
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as  the  constant  value  of  range  between  the  point  source  and  the  discrete  point  scatterer  when  the 
transmitted  acoustic  field  is  first  incident  upon  the  discrete  point  scatterer  at  time  instant  t'  after 
motion  begins  at  time  instant  tm  where  t'  >  tm.  The  range  r01  is  given  by  (2.1-26). 


Now  let  us  solve  for  the  constant  value  of  range  |r"2|.  Substituting  (2.3-153)  and  (2.3-159) 

into  (2.3-61);  (2.3-153),  (2.3-154),  and  (2.3-159)  into  (2.3-62);  (2.3-161)  into  (2.3-63);  (2.3-162) 
into  (2.3-64);  (2.3-153),  (2.3-154),  (2.3-159),  and  (2.3-161)  into  (2.3-65);  and  (2.3-166)  into  (2.3- 
66)  yields 


A\  —  2.  rl2  (^i,2#  V0), 

(2.3-167) 

a2=v0\ 

(2.3-168) 

A3  =  2ri,2(«i,2*Vo). 

(2.3-169) 

■t* 

II 

(2.3-170) 

A5  =  2V02, 

(2.3-171) 

(2.3-172) 

With  the  use  of  (2.3-57)  and  substituting  (2.3-170)  into  (2.3-58);  (2.3-169),  (2.3-171),  and  (2.3- 
172)  into  (2.3-59);  and  (2.3-167),  (2.3-168),  and  (2.3-172)  into  (2.3-60),  we  obtain 


L"  _^o±V 

'^o+4^o^o 

1 r  1,2  - 

2^o 

where 


(2.3-173) 


and 


v 2 

=  l--2- 
~0  1  2  » 
C 


So=2rfA^A+2r, 


0.1 


(2.3-174) 

(2.3-175) 


^=r^  +  2r  r  Cgy* Yp)  +  2 

(so  —  r0,l  2  ^zr0,1ri,2  +  '1,2 ' 

c  c 


(2.3-176) 


The  ranges  r0>1  and  rli3  are  given  by  (2.1-26)  and  (2.1-27).  respectively.  The  solution  given  by  (2.3- 
173)  is  the  constant  value  of  range  between  the  discrete  point  scatterer  and  the  receiver  when  the 
scattered  acoustic  field  is  first  incident  upon  the  receiver  at  time  instant  t  after  motion  begins  at  time 
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instant  tm  where  t>t'>tm.  The  decision  to  use  either  the  plus  or  minus  sign  in  (2.3-173)  is 
dictated  by  the  fact  that  range  must  be  positive. 

Let  us  solve  for  the  time-varying  range  r "(/)  next.  Substituting  (2.3-168),  (2.3-170),  and 

(2.3-171)  into  (2.3-72);  (2.3-167)  through  (2.3-169)  and  (2.3-171)  into  (2.3-73);  and  (2.3-167) 
and  (2.3-168)  into  (2.3-74)  yields 

/*  =  1,  (2.3-177) 

.ff(0  =  0,  (2.3-178) 

and 

m  =  Vj (At)2  +  2r12(n12.  \0)At  +  r?2,  (2.3-179) 

where  [see  (2.3-75)  and  (2.3-76)] 

A  t>tm  +  r,  (2.3-180) 

and 


(2.3-181) 

is  the  time  delay  in  seconds  (the  amount  of  time  it  takes  for  the  transmitted  acoustic  signal  to  begin 
to  appear  at  the  receiver  after  motion  begins  at  time  instant  tm )  where  the  constant  values  of  range 

|ro.i|  and  |r,':2|  are  given  by  (2.3-166)  and  (2.3-173),  respectively.  Substituting  (2.3-177)  through 
(2.3-179)  into  (2.3-71)  yields 

(2.3-182) 

Now  let  us  solve  for  the  time-varying  range  |rQ>t(0|.  With  the  use  of  (2.3-78),  we  can  write 
that 


(2.3-183) 

where  At  is  given  by  (2.3-180),  |r"(r)|  is  given  by  (2.3-182),  and  r  is  given  by  (2.3-181). 

Let  us  solve  for  the  time-varying  unit  vector  n'Q  l(t)  next.  If  we  substitute  (2.3-159)  into  (2.3- 
82),  then 
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«Sj(0= 


(2.3-184) 


t>tm  +  T. 

m 


Equation  (2.3-184)  does  not  make  physical  sense  because  we  need  the  time-varying  vector  r^,(r)  in 
the  numerator  and  not  the  constant  vector  r0 , .  Therefore,  we  need  to  derive  an  equation  for  r'0 u(f) . 
Before  we  begin  the  derivation,  let  us  note  that  although  (2.3-184)  does  not  make  physical  sense,  if 
we  substitute  (2.3-159)  into  (2.3-17),  then 


ro,i —  ro,i 


(2.3-185) 


and,  as  a  result. 


l0,l 


'0,1- 


(2.3-186) 


Equation  (2.3-186)  does  make  physical  sense  and  it  agrees  with  (2.3-166). 

We  begin  the  derivation  of  the  equation  for  the  time-varying  vector  r'  ,(r)  by  substituting 
(2. 1  -6)  into  (2.3- 1 85)  yielding 


ro.,=  r>-r0- 


(2.3-187) 


We  then  generalize  (2.3-187)  as  follows: 

=  t>tm  +  r,  (2.3-188) 


where  1Ra(t)  is  the  time- varying  position  vector  from  the  origin  to  the  transmitter  given  by  (2.3-4), 
and  J£j(f)  is  the  time- varying  position  vector  from  the  origin  to  the  discrete  point  scatterer  given  by 
(2.3-5).  Substituting  (2.3-4),  (2.3-5),  and  (2.3-158)  into  (2.3-188)  yields 


ro,i(0  =  ro,-lt~(L  +  Vo,  t>tm  +  z. 


(2.3-189) 


where  r01  is  given  by  (2.1-6)  and  t  is  given  by  (2.3-181).  Note  that  if  we  evaluate  (2.3-189)  at 
t=tm  + 1,  then 

m  7 


•*0,1  T)  —  r0,l  > 

where  r'  ,  is  given  by  (2.3-185).  Therefore, 


(2.3-190) 


(2.3-191) 


where  |r^  ,|  is  given  by  either  (2.3-186)  or  (2.3-166).  Equation  (2.3-191)  agrees  with  the  general 
result  given  by  (2.3-80). 

As  a  result  of  the  above  analysis,  the  time-varying  unit  vector  n'a  ,(t)  is  given  by 
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(2.3-192) 


t>tm  +  T, 


where  i*o,(f)  is  given  by  (2.3-189)  and  r  is  given  by  (2.3-181).  And  by  referring  to  (2.3-85)  and 
(2.3-86),  we  can  write  that 


and 


qjj(f)  =  cos"' wjj(0,  t>tm  +  T, 


(2.3-193) 


f 

i/^o.i(0  =  tan"1 

V 


t  >  /  +  Z, 


(2.3-194) 


where  u'0l(t),  v'0  l(t),  and  Wg/t)  are  the  dimensionless,  time-varying  direction  cosines  with  respect 
to  the  X,  Y,  and  Z  axes,  respectively,  associated  with  the  time-varying  unit  vector  n'  ,(t)  given  by 
(2.3-192),  and  r  is  given  by  (2.3-181).  Equations  (2.3-193)  and  (2.3-194)  are  the  angles  of 
incidence  at  the  discrete  point  scatterer  that  are  to  be  used  to  evaluate  die  scattering  function  of  the 
discrete  point  scatterer  given  by  (2.3- 1 0 1 ). 

Let  us  solve  for  the  time- varying  unit  vector  n"2(t)  next.  If  we  substitute  (2.3-153),  (2.3-159), 
and  (2.3-161)  into  (2.3-88),  then 


where 


r"(t) 

KiW^rr— i’ 

t  >  tm  +  T, 

|«u«)| 

(2.3-195) 


(2.3-196) 


r12  is  given  by  (2.1-14),  At  is  given  by  (2.3-180),  and  z  is  given  by  (2.3-181).  Note  that  (2.3-196) 

can  be  obtained  by  substituting  (2.3-13),  (2.3-25),  (2.3-75),  (2.3-153),  (2.3-159),  and  (2.3-161)  into 
(2.3-28).  And  by  referring  to  (2.3-91)  and  (2.3-92),  we  can  write  that 


(2.3-197) 


and 
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(2.3-198) 


where  u"2(t),  v"2(f) ,  and  w"2{t)  are  the  dimensionless,  time-varying  direction  cosines  with  respect 
to  the  X,  Y ,  and  Z  axes,  respectively,  associated  with  the  time- varying  unit  vector  n"2(t)  given  by 
(2.3-195),  and  t  is  given  by  (2.3-181).  Equations  (2.3-197)  and  (2.3-198)  are  the  angles  of  scatter 
at  the  receiver  that  are  to  be  used  to  evaluate  the  scattering  function  of  the  discrete  point  scatterer 
given  by  (2.3-101). 

Finally,  let  us  solve  for  the  time- varying  position  vector  r 2{t)  from  the  origin  to  the  receiver. 
If  we  substitute  (2.3-153)  and  (2.3-161)  into  (2.3-95),  then 


r"(0  =  r2  +  ArV0,  t>tm  +  T, 


(2.3-199) 


where  A t  is  given  by  (2.3-180)  and  t  is  given  by  (2.3-181). 

In  order  to  fully  appreciate  the  importance  of  the  results  contained  in  this  S  AS  example,  keep 
in  mind  that  they  are  based  on  the  exact  time-varying  ranges  between  the  transmitter  and  discrete 
point  scatterer,  and  between  the  discrete  point  scatterer  and  receiver,  and  2)  the  exact  time-varying 
angles  of  incidence  at  the  discrete  point  scatterer,  and  the  exact  time-varying  angles  of  scatter  at  the 
receiver  that  are  used  to  evaluate  tire  frequency  dependent  scattering  function  of  the  discrete  point 
scatterer.  In  contrast,  Bonnifant  [8],  for  example,  1)  initially  ignores  the  height  of  the 
transmit/receive  platform  above  the  ocean  bottom  when  computing  ranges,  2)  he  assumes  that  the 
target’s  reflectivity  is  constant  -  not  a  function  of  frequency  and  angles,  and  3)  that  the 
transmit/receive  platform  is  stationary  during  signal  transmission  and  reception  -  the  common 
“stop  and  hop”  assumption  that  he  later  attempts  to  correct  for  with  an  approximate  phase  factor 
correction.  We  have  made  no  such  assumptions  in  this  example  and,  as  a  result,  no  corrections  are 
necessary. 
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3  Summary 


The  complex  frequency  response  of  the  ocean  was  derived  for  the  following  three  different 
bistatic  scattering  problems:  1)  no  motion,  2)  only  the  discrete  point  scatterer  is  in  motion,  and  3)  all 
three  platforms  (the  transmitter,  discrete  point  scatterer,  and  receiver)  are  in  motion.  The  propagation 
of  the  small-amplitude  acoustic  signals  in  the  ocean  involved  in  the  bistatic  scattering  problems  was 
treated  as  transmission  through  a  linear,  time-variant,  space-variant  filter. 

Scatter  from  a  discrete  point  scatterer  was  modeled  via  the  scattering  function,  which  is  a 
complex  function  (magnitude  and  phase)  and  is,  in  general,  a  function  of  frequency,  the  direction  of 
wave  propagation  from  the  source  to  the  scatterer,  and  the  direction  of  wave  propagation  from  the 
scatterer  to  the  receiver.  In  addition  to  the  scattering  function,  frequency-dependent  attenuation  was 
taken  into  account  in  order  to  model  the  propagation  of  sound  from  transmitter  to  discrete  point 
scatterer,  and  from  discrete  point  scatterer  to  receiver. 

The  speed  of  sound  and  ambient  density  of  the  ocean  were  treated  as  constants.  Therefore, 
sound  rays  will  travel  in  straight  lines.  We  only  concerned  ourselves  with  solving  for  the  direct  ray 
path  between  transmitter  and  discrete  point  scatterer,  and  from  discrete  point  scatterer  to  receiver.  As 
a  result,  the  three  platforms  were  treated  as  being  in  an  unbounded,  homogeneous  ocean  medium. 
However,  sound  propagation  between  the  transmitter  and  the  ocean  surface  and  bottom,  and  from 
the  ocean  surface  and  bottom  to  the  receiver  can  be  handled  in  the  same  way  as  was  done  for  the 
discrete  point  scatterer. 

Section  2.1  was  devoted  to  the  first  bistatic  scattering  problem,  which  involves  no  motion  -  the 
transmitter,  discrete  point  scatterer,  and  receiver  are  not  in  motion.  The  exact  solutions  for  the 
angles  of  incidence  at  the  discrete  point  scatterer  and  the  angles  of  scatter  at  the  receiver  were 
derived.  An  example  was  worked  out  at  the  end  of  Section  2.1  showing  how  the  general  bistatic 
scattering  results  reduced  for  a  monostatic  ( backscatter )  scattering  geometry  for  the  no  motion 
case. 

Section  2.2  was  devoted  to  the  second  bistatic  scattering  problem  when  only  the  discrete  point 
scatterer  is  in  motion.  Motion  was  allowed  to  start  at  an  arbitrary  time  instant  tm  seconds  as 
opposed  to  zero  seconds.  Two  new  major  results  were  presented  in  Section  2.2:  1)  the  exact  time- 
varying  ranges  between  the  transmitter  and  discrete  point  scatterer,  and  between  the  discrete  point 
scatterer  and  receiver  were  derived,  and  2)  the  exact  time-varying  angles  of  incidence  at  the  discrete 
point  scatterer,  and  the  exact  time-varying  angles  of  scatter  at  the  receiver  were  also  derived.  An 
example  was  worked  out  at  the  end  of  Section  2.2  showing  how  the  general  bistatic  scattering 
results  reduced  for  a  monostatic  ( backscatter )  scattering  geometry  for  the  case  when  only  the 
discrete  point  scatterer  is  in  motion. 

Section  2.3  was  devoted  to  the  third  bistatic  scattering  problem  when  all  three  platforms  are  in 
motion.  Motion  was  allowed  to  start  at  an  arbitrary  time  instant  tm  seconds  as  opposed  to  zero 
seconds.  Two  new  major  results  were  presented  in  Section  2.3:  1)  the  exact  time-varying  ranges 
between  the  transmitter  and  discrete  point  scatterer,  and  between  the  discrete  point  scatterer  and 
receiver  were  derived,  and  2)  the  exact  time-varying  angles  of  incidence  at  the  discrete  point 
scatterer,  and  the  exact  time-varying  angles  of  scatter  at  the  receiver  were  also  derived.  Three 
examples  were  worked  out  at  the  end  of  Section  2.3.  The  first  example  showed  that  the  exact  results 
derived  in  Section  2.3  reduced  to  the  exact  results  derived  in  Sections  2.1  and  2.2  when  appropriate 
values  were  used  for  the  various  parameters.  This  is  a  very  important  example  because  it  validates 
the  correctness  of  the  general  solution  derived  in  Section  2.3.  The  second  example  showed  how  the 
general  bistatic  scattering  results  reduced  for  a  monostatic  {backscatter)  scattering  geometry  for  the 
case  when  all  three  platforms  are  in  motion.  The  third  example  showed  how  die  general  bistatic 
scattering  results  can  be  applied  to  a  synthetic  aperture  sonar  (SAS)  system  trying  to  image  a 
nonmoving  target  on  the  ocean  bottom  without  having  to  make  several  common  simplifying 
assumptions. 

Finally,  it  is  important  to  note  that  for  problems  involving  motion,  the  solutions  for  the  exact 
time-varying  ranges  between  the  transmitter  and  discrete  point  scatterer,  and  between  the  discrete 
point  scatterer  and  receiver  derived  in  this  report  are  also  valid  in  an  inhomogeneous  ocean  where 
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the  speed  of  sound  and  ambient  density  are  functions  of  position  since  solving  for  a  range 
represents  a  problem  in  mechanics  not  wave  propagation.  However,  travel  times  and  angles  of 
incidence  and  scatter  are  different  in  an  inhomogeneous  ocean  compared  to  a  homogeneous  ocean 
because  of  the  complicated  trajectories  of  sound  rays  in  an  inhomogeneous  ocean. 
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